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Photoelectric  data  in  the  UDV  systeei  were  obtained  by  the  author 
in  1975  and  1975  at  Bosemary  Hill  Observatory,  These  data  of  RS  Canum 
Venatlconen  were  analyzed  together  with  the  published  Catania  data.  The 
analysis  Included  the  use  of  a modified  version  of  0.  8.  Hood's  (A.J.  ]i, 
701]  computer  program  HINA.  A truncated  Fourier  series  representation 
of  the  distortion  wave  was  subtracted  from  the  observed  light  curve  so 
improved  elements  could  be  determined  from  Che  difference  by  use  of  the 
KINK  program.  This  process  was  not  at  all  like  the  usual  rectification 
procedure;  it  was  an  iterative  procedure  which  greatly  improved  the 
confidence  in  the  determination  of  tne  elements  of  RS  CVn.  Hith  the 
single  exception  of  the  effective  tHnperature  of  the  KO  III  star,  the 
eclipsing  elements  are  consistent  for  all  epochs  and  colors  studied. 

Using  these  elements,  it  was  pcssiple  for  the  distortion  wave  to  be 
isoleted;  therefore,  it  will  now  be  possible  for  theories  of  its  origin 
and  structure  to  be  tested  more  critically.  The  results  ere  consistent 


vritti  a spotted  surface  for  tfle  KD  star.  Hcwever,  tlie  sunspot  analogy 
does  not  accurately  describe  tbe  nature  of  tbe  spots. 


CHAPTER  I 
INTROOUCnW 

Topic  of  Current  Resaarch 

The  eclipsing  binary  star  RS  Canuoi  Venatlconan  (RS  CVn)  1s  the 
prototype  of  a very  interesting,  but  poorly  understood,  group  of  stars. 

It  nas  Deen  chosen  for  further  study  because  It  has  the  most  complete 
observational  record  of  this  group,  and  because  It  Is  sufficiently 
bright  to  allow  observations  with  the  available  equipment.  The  RS  CVn 
binaries  [also  known  In  the  literature  as  RS  CVn  systems,  RS  CVn-type 
binaries,  RS  CVn  stars,  etc.)  nave  been  recently  reviewed  by  Hell  [1R76}. 
Ke  concludes  that  the  group  1s  best  defined  by  systems  in  which  (1)  the 
hbtler  star  is  of  spectral  type  P or  G and  luminosity  type  V or  IV, 

(2]  strong  Ca  II  emission  Is  present  In  the  H and  K lines,  and  (3)  the 
orbital  period  of  the  binary  system  1s  between  approximately  a day  and 
2 weeks.  There  are  other  characteristics  which  are  comoon  to  many  of  the 
systems  listed  by  Hall,  but  which  were  deemed  unnecessary  for  the  purpose 
of  classification.  Some  of  these  characteristics  would  be  mass  ratios 
near  unity,  cooler  components  near  spectral  class  KD  III,  and  photo- 
metric distortions  superimposed  on  the  normal  eclipse  phenomena  (if  It 
exists).  It  Is  the  complications  caused  ty  the  last  of  these  characteris- 
tics to  which  this  dissertation  1s  addressed. 

A new  set  of  observations  of  RS  CVn  will  be  presented  which 
will  extend  the  basic  astronomical  khowledge  of  the  system  and  confirm 
the  earlier  observations.  A new  technique  of  solution  Is  developed 


which  drcanvents  the  previous  problems  encountered  when  derivinj  ele- 
ments free!  a light  curve  with  unejplained  light  variations.  A natural 
result  of  this  technique  is  to  isolate  the  unexplained  light  variation 
(referred  to  as  a 'distortion  wave'  in  this  work)  so  that  It  may  be 
studied  Independently,  It  Is  also  hoped  that  the  times  of  minima  can 
be  improved  with  this  technique  so  that  a better  represehtatlon  of  the 
period  viriations  can  be  obtained. 

In  the  remainder  of  this  chapter  the  literature  on  RS  CVn  and 
RS  CVn  binaries  will  be  reviewed  and  the  scope  of  the  dissertation  will 
be  outlined.  In  the  chapters  to  follow  the  new  photoelectric  data  will 
be  presented,  starting  with  the  instrumentation  and  procedure  used. 
Following  this  will  be  a discussion  of  the  distortion  wave  and  period 
variations.  The  dissertation  will  conclude  with  a summary  of  the 
results,  interpretation  as  to  their  meaning,  and  suggestions  for  future 
research. 

Historical  Review 

The  review  of  the  litereture  shall  start  with  a chronological 
history  of  the  research  on  RS  CVn.  Khen  this  review  reaches  the  early 

RS  Canu  Venaticorum  from  1914  to  19?2 

RS  CVn  was  first  reported  as  a variable  star  by  L.  Ceraski  in 
1914.  She  reported  that  the  star  was  observed  to  be  fainter,  on 
photographic  plates  read  by  M.S.  Blaaka.  on  7 April  1896  (by  1?5), 

IS  May  1907  (by  l^O).  and  3 May  1911  Cby  1^1.  Numerous  observers 
studied  this  systen  during  the  following  16  years.  The  results  of 


these  observetions  culmlnetee  In  the  solution  of  the  light  curve  presented 
by  Sitterly  [1630].  Sitterly  surmarized  the  considerable  observa> 
tional  naterial  obtained  between  1914  and  1923  by  Townley  (19IS), 
Hoffmeister  (IdlSi  1919]i  Haggini  (1916],  Gadcinskl  (1926)i  Schneller 
(1923)  and  others.  He  presented  the  tie«s  of  isinima  that  these  various 
observers  had  collected  up  to  1923  and  discussed  the  period  changes. 

He  noted  that  the  period  steadily  Increased  during  the  interval  1900 
to  1928-  This  can  easily  be  verified  from  the  times  of  minina  as  plotted 
on  the  0*0  diagram  [Observed  time  of  minima  minus  the  Calculated  time  of 
minima  from  a linear  ephedieris].  He  had  accumulated  710  visual  estimates 
using  the  23-inch  Princeton  refractor  and  351  photographic  astimatas  from 
Harvard  plates.  The  visual  estimates  were  obtained  during  the  years  1920 
to  1922  [mostly  in  1921).  while  the  photographic  estimates  covered  a much 
longer  time  (from  1892  to  1922).  The  visual  light  curve  showed  asynmetry 
outside  eclipse  and  in  primary  eclipse.  The  photographic  light  curve,  on 
the  other  hand,  showed  no  asynmetry  outside  eclipse,  nor  did  it  show  a 
secondary  eclipse.  In  retrospect,  this  can  be  understood  in  terms  of  a 
moving  distortion  wave  averaging  out  the  asynmetry  over  the  much  longer 
Interval  of  time  during  which  the  photographic  estimates  were  oiace.  This 
process  would,  of  course  result  in  a much  larger  standard  error  for  points 
of  the  photographic  light  curve,  but  this  would  be  expected  In  any  case 
because  of  Che  lower  precision  of  this  technigue.  Except  for  the  asymme- 
try outside  eclipse,  Sitter1y‘s  visual  light  curve  solution  was  in  good 
agreement  with  the  photographic  light  curve.  The  solution  was  CMUJlicated 
by  a small  secondary  aclipae,  only  0?03  deep,  which  ccmpares  with  a depth 
of  o’i'2  by  modern  photoelectric  photometry.  Sitterly  decemlned  a depth  of 


for  primary  eclipse,  which  is  in  reasonable  agresiiBPt  with  a modern 
value  of  iTl.  he  attempted  a solution  in  which  the  following  sides  of 
both  stars  were  more  luminous  chan  the  advancing  sides.  This  was  neces- 
sary in  order  to  eaplain  the  outside  eclipse  asynnetry,  which  amounted 
to  a peak  to  peak  difference  of  0?13  between  the  portion  of  the  light 
curve  preceding  and  following  primary  eclipse.  This  gave  much  better 
residuals  than  the  other  solutions  but  it  was  discarded  because  it 
seemed  so  physically  inekplicable.  particularly  since  It  required  the 
fainter  star  to  have  a difference  in  brightness  from  one  side  to  the 
other  of  55S1  The  solution  which  was  finally  adopted  as  the  best 
representation  of  the  true  values  was  obtained  by  solving  the  visual 
light  curve  without  any  rectification.  The  resulting  elements  were: 


Uniform  disks  assumed. 


where  the  subscript  b (or  2)  refers  to  the  brighter  star,  and  the  sub- 
script f (or  1}  refers  to  the  fainter  star.  The  estimated  probable 
error  for  a visual  observation  was  f o'i’oasa,  while  the  estimated 
probable  error  for  a photographic  observation  was  + 0?Q9.  It  is  there- 
fore wise  that  he  chose  to  base  the  solution  on  the  visual  light  curve. 
The  comparison  star  used  was  B&r36°  2345,  and  the  differential  magni- 
tudes were  originally  given  as  comparison  minus  variable. 

In  the  same  year,  Joy  (1930)  published  spectroscopic  data  on 
R5  CVn.  He  used  the  results  of  Sitterly’s  photometric  solution  along 
with  his  own  spectroscopically  determined  radial  velocities  to  arrive 


at  aOtoluU  dimensions.  A sumnary  of  Joy's  results  Is  given  in  Table  1 
(the  subscripts  have  the  same  meaning  as  noted  above). 


JOY'S  ABSOLUTE  DIKMIQNS 


ICj.91.6  km/s 

rijSin^l  ■ 1,79 

kj  ■ 99.0  kla/s 

mjSin^  = 1,66 

V • -8,9km/s 

mj/m,  = 0.93 

fb  ■ “o 

.jsin  i . 6,53  X loS  km 

'■f  • 5.3  l!„ 

m^  w 1.85 

• 1.71 

Joy  concluded  from  the  spectroscopic  absolute  magnitude  and  from 
the  spectral  class  that  both  stars  are  dwarfs.  However,  the  siae  and 
brightness  of  the  fainter  star  were  in  direct  conflict  with  this.  Joy 
ended  his  article  by  saying,  ".  . , the  secondary  star  is  certainly  pecul- 
iar and  cannot  readily  be  classified  among  other  stars  for  which  we  know 
the  physical  characteristics"  (p.  45).  Joy  agreed  with  Sitterly's  inter- 
pretation that  the  fainter  secondary  star  is  "a  sort  of  'subgiant.'" 

C.  Payne-Gaposchkin  (ig39)  made  a study  of  the  period  changes  and 
light  curve  asymoetries  Of  AS  CVn  as  a part  of  a paper  on  variable  stars. 
She  collected  approximately  4000  photographic  estimates  of  the  brightness 
from  Hsrvard  plates  exposed  between  the  years  1892  and  1938.  From  the 
data,  17  rimes  of  minima  were  computed  and  plotted.  In  addition,  a light 
curve  was  made  by  combining  individual  points  Into  normals.  In  primary 
eclipse  each  normal  was  mede  up  of  10  individual  points,  while  outside 
eclipse  each  nonnal  was  cpmprised  of  IQO  individual  points.  The 
nomml  points  were  computed  after  each  individual  point 


was  correctad  for  the  change  in  period,  as  she  had  estimated  it.  The 
resulting  light  curve  agreed  with  Sicterly's  in  its  general  characteris- 
tics. However,  it  did  not  agree  with  the  asjnmetry  that  Sitterlj  had 
observed  from  his  visual  light  curve.  Recall  that  Sitterly  had  pre- 
sented a visual  and  photographic  light  curve  for  115  CVn,  but  his  photo- 
graphic light  curve  failed  tp  show  a secondary  minimum  or  any  asymmetry 
outside  eclipse.  Payne-Gaposchkin,  pn  the  other  hand,  found  a shallow 
secondary  eclipse  in  her  photographic  light  curve  (pf  depth  hl'06},  and 
a difference  in  Che  height  between  Che  two  navima  of  only  a few  hun- 
dredths of  a magnitude.  She  noted,  however,  that  when  the  light  curve 

scmewhat  larger  difference  was  found.  In  an  effort  to  summarize  the 
differences  in  magnitude  between  the  two  maxima,  she  presented  a table 
pf  these  differences  at  four  epochs  of  minima-  She  noted  that  if  the 
difference  is  assumed  to  Pe  a periodic  phenomena,  then  the  period  would 
be  approximately  BOOO  days  (22  years)  or  1670  orbital  periods,  i.e-.  the 
approximate  period  of  the  orbital  period  changes.  She  estimated  that 
an  absolute  orbital  radius  of  2 X 10^  Ion  would  be  necessary  to  explain 
the  period  changes  as  a llght-time  effect.  Further,  she  determined 
that  if  this  were  due  to  a third  body,  then  the  required  mass  of  4 fi 
would  seem  larger  than  could  possibly  go  undetected  photcmetrically. 

This  represented  the  first  convincing  evidence  that  Che  period  changes 
were  of  an  unusual  nature,  and  the  first  time  that  the  differences  in 
heights  of  maxima  were  ascribed  to  be  of  a possibly  periodic  nature. 

The  most  significant  aspect  of  this  paper  wes  that  the  differences  in 
the  height  of  the  two  raxiraa  were  ascribed  to  the  redder  (and  fainter) 


of  tn«  two  stars.  This  was  basad  on  her  realization  that  the  size  of  the 
asyinnetry  was  greater  1n  the  visual  light  curve  than  1n  the  blue  photo- 
graphic light  curve. 

H.  A.  KlUner  (19A7)  noted  that  a define  group  of  stars  eiisted 
which  showed  Ca  II  emission.  In  this  short  note  he  collected  the 
scattered  data  on  the  13  systesB  known  to  him  at  that  tine  (RA  Eri.  SS  Can, 
AR  hon,  RU  Cnc,  RS  Cnc,  RW  UHa,  RS  CVn,  SS  Boo.  UW  Dra,  Z Her.  M Her, 

RT  Lac  and  AR  Lac).  Hlltner  found  from  his  study  of  RS  CVn  the  same 
characteristics  reported  by  other  observers  whose  data  he  had  collected 
together.  The  emission  was  noted  to  be  usually  from  the  fainter  of 
the  two  stars  1n  the  eclipsing  system  and  the  emission  was  apparently 
".  . . not  situated  uniformly  around  the  parent-star"  {p.ABI).  In  fact 
he  stated  that  the  emission  approached  invisibility  only  at  the  secondary 
eclipse  when  the  fainter  star  Is  partially  covered  by  the  brighter,  but 
smaller  star.  This  led  h1n  to  agree  with  Struve  (1544)  and  others  thel 
the  Ca  II  emission  In  these  stars  originates  from  a tidally  extended  re- 
gion above  the  photosphere  of  the  fainter  star.  The  emitting  region  was 
isolated  at  the  extremes  of  this  tidal  bulge  along  the  line  Joining  the 
two  stars.  Struve's  and  Hlltner's  papers  could  be  considered  Che  first 
recognition  that  there  was  a group  of  stars  (today  referred  to  as 
RS  CVn  binaries)  which  warranted  consideration  as  a group. 

SraCton  (19SO)  studied  19  systens  which  showed  H and  < lines  of 
Ca  I!  emission  In  the  last  section  of  a paper  published  just  a few  years 
after  Hlltner's.  This  list  Included  the  13  stars  In  Hlltner's  paper 
plus  six  others  which  were  not  eclipsing  binaries,  but  were  of  a binary 
six  stars  were  primarily  more 


luminous  and  of  longer  periods 


tlic^e  in  Kiltner'&  list.  Gratton  attempted  to  give  additional 
evidence  that  the  tidal  beige  proposed  by  Struve  and  fiUtner  Mas  the 

Briefly,  Grattan's  analysis  folloMed  these  lines.  By  assuming 
that  all  stars  tvith  a tidal  distortion  above  some  lower  limit  would 
show  emission  lines,  it  was  then  expected  that  the  giant  stars  which 
show  emission  would  have  longer  orbital  periods.  This  condition  led 
to  the  relation 


where  it  was  assumed  that  the  variation  in  mass  could  be  neglected,  R 
was  the  radius  of  the  star,  a^  and  ^2  the  relative  orbital  radii, 
and  was  a constant.  Since  the  stans  being  studied  were  of  approxi- 
mately the  same  spectral  type  (and  therefore  had  the  same  surface  bright- 
ness), the  radius  was  proportional  to  the  square  root  of  the  luminosity. 
In  addition,  the  sum  of  the  relative  orbital  radii  was  proportional  to 
the  period  (P)  to  the  two-thirds  power.  Therefore  the  period  and  the 
absolute  magnitude  CH)  for  the  binaries  showing  H and  K emission  lines 
were  above  a line  0.3  n r log  P * constant. 

Gratton  plotted  the  absolute  magnitude  versus  log  P for  13  of 

AR  lac,  and  RO  Cue  were  excluded  because  of  the  lack  Of  data.  Gratton 
stated,  "It  is  seen  that  there  is  a definite  correlation  between  M and 
P;  wo  may  tentatively  take  the  limiting  line  0.3  M + log  P “ 2.5"(p.  40). 
Stars  above  the  line  0.3  fi  t log  P ■ 1.0  are  practically  in  physical 
contact;  for  these  stars  the  distortion  would  be  too  great  and  might 
cause  instability. 


present  writer 


tnat  the  assumptions  wfiich  led  to 
this  conclusion  are  reasonable,  with  the  possible  exception  of  the  cor- 
relation between  the  eml5SlOh  ahd  the  degree  of  tidal  distortion.  Here 
recently,  Young  and  Konlges  (1977)  have  Investigated  this  relationship; 
this  will  be  discussed  later, 

Keller  and  Limber  (1951)  used  an  unfiltered  1P21  ohotomultlpller 
tube  to  obtain  a light  curve  of  RS  CVn  during  the  spring  of  1949.  They 
used  60*35°  2421.  60*35°  2422.  and  Bl>*35°  2418  as  comparison  stars,  with 
the  differential  lisninosity  expressed  relative  to  60*35°  2422.  Inside 
eclipse.  Individual  points  were  tabulated,  but  outside  eclipse  only  the 
normal  points  were  given.  Frpm  an  expanded  plot  of  these  data  around 
the  priiiary  minimum  Keller  and  Limber  noticed  a systenatic  difference  in 
luminosity.  They  found  that  this  variation  in  the  depth  of  primary  was 
most  probably  due  to  Intrinsic  variability  of  tbe  fainter  star.  They 
also  noted  the  same  asynoetry  outside  eclipse  observed  by  SHterly  and 
Payna-Gaposchkin.  It  Is  Important  to  note  that  the  observations  of 
Keller  and  Limber  were  taken  over  only  a three-month  Interval,  so  that 
the  nomels  they  produced  outside  eclipse  were  not  affected  as  much  by 
variations  in  the  asymmetry  as  was  the  case  for  the  observations  of 
SUterly  (mostly  one  year  of  data)  or  Payne-Saposchkln  (nearly  47  years 
of  data).  The  asyimetry  observed  by  Keller  end  Limber  was  of  the  same 
general  character  as  that  of  the  previous  observers;  It  had  a higher 
nHxImum  after  primary  eclipse  then  prior  to  primary  eclipse.  The  dif- 
ference In  the  two  maxima  was  0*!’04.  The  rectified  light  curve  was 
solved  using  a Heib  darkening  coefficient  of  0.8  with  the  aid  of  the 
Herrill  (1960)  tables.  The  elanents  so  detemlned  «re; 


• 0.726,  Lf  ■ 0-274, 


frotn  the  lunlnosHy  and  radius  of  each  component  they  computed 
the  ratio  of  surface  intensities.  This  ratio,  aporoximately  20,  could 
not  be  reconciled  eith  Che  value  of  5.5  they  conpuced  on  Che  basis  of 
the  eneroy  distribution  of  the  stars  Ccorrected  for  the  effects  of  the 
Balmer  discontinuity,  atmospheric  transmission,  and  the  response  of 
the  instnmentation).  This,  of  course,  was  the  same  problem  Joy  had 
oomnented  upon. 

Keller  and  Limber  very  crudely  estimated  the  poiytropic  inOea 
of  Che  fainter  star  to  be  2.65.  They  used  an  0-C  diagram  made  from  all 
the  available  times  of  minima,  and  an  orbital  eccentricity  determined 
from  a single  estimate  of  the  time  of  secondary  minimum  based  on  their 


While  it  is  obvious  that  their  ohservations  are  the  best 
available  at  that  tine.  It  is  equally  obvious  Cas  can  be  seen  fnc 
light  curve  in  their  oaper)  that  the  time  of  secondary  minimum  is 
greatly  complicated  by  the  asyametry  of  the  light  curve.  It  ms 
fore  wise  of  them  to  state,  "It  is  not  possible  to  obtain  the  ohe 
the  secondary  with  much  precision'.'  (p.  647). 

Bidelman  C16S4)  published  a list  of  stars  which  were  knov 
show  emission  lines  and  whose  spectral  types  were  later  than  B. 
list  for  stars  showing  Ca  il  emission,  «h 
binanes,  continues  to  be  referenced  coda; 
the  other  stars  in  Hiltner's  paper  were  ai 


spectroscopic 
s well  es  all  of 
stars  in  Bidelman's 


Popper  (1961}  observed  RS  CVn  lo  an  effort  to  Improve  the 
quality  Of  the  spectroscopic  orbit  and  to  add  photoelectric  observations 
on  a standard  systetn  to  the  body  of  existing  data.  He  oOtained  12  spec- 
trograms between  the  years  1954  and  1958  which  showed  the  same  general 
features  previously  observed  by  Joy  (1930)  and  Miitner  (194?).  The  only 
major  difference  between  his  spectrograms  and  those  of  the  previous 
observers  was  the  sharpness  of  the  absorption  lines  due  to  the  brighter 
component  of  spectral  class  F4.  Joy  had  designated  the  spectral  type  as 
F4n,  which  would  mean  that  the  lines  were  "nebulous."  Popper  noted  the 
bright  H and  X emission  lines  of  Ca  II  from  the  cooler  (fainter)  compo- 
nent, as  well  as  the  variable  emission  intensity  of  the  H line.  The 
solution  of  Popper's  radial  velocity  curve  was  compared  with  Joy's. 

6otn  observers  assumed  a circular  orbit.  Popper's  results  gave  slightly 
smaller  masses  for  both  stars.  Ih  addition,  he  believed  that  the  differ- 
ence in  the  systemic  radial  velocity  of  6 Icm/s  was  probably  a real 
(but  unexplained)  phenomenon.  The  photcnetry  was  done  during  the  suncner 
of  1956,  in  January  1958,  and  in  January  1959.  No  attenpt  was  made  to 
cover  the  entire  light  curve.  The  ccmparison  star  was  KD  1I477S 
(BO»35°  2420),  i.e.,  the  same  Star  used  later  by  Chisari  and  Ucona 
(1965),  Catalano  and  Rodono  (1967)  and  the  present  writer.  The  photom- 
etry of  this  star  from  the  three  seasons  resulted  in  the  following 
magnitudes: 

V = 8?42,  B-V  • ♦0?46,  and  U-B  = -0%4. 

This  star  was  0?03  fainter  in  1956  then  it  was  during  the  1956-59 
Observing  seasons.  Popper  was  npt  Sure  whether  this  difference  was  due 
to  a real  variation,  to  the  difference  in  the  zenith  angles  involved, 
or  to  the  difference  in  the  set  of  standard  stars  used. 


As  expected,  primary  eclipse  was  oCserved  to  be  deeper  for 
sborter  wavelengths,  and  secondary  mlnimun  was  observed  to  be  deeper  In 
the  yellow  than  it  was  In  the  blue.  Surprisingly,  the  depth  of  prieary 
eclipse  was  o'I's  different  in  19S6  than  it  was  in  19S8-S9.  This  must  be 
attributed  to  the  fainter  star,  since  this  eclipse  is  a total  one.  ho 
such  variation  was  found  in  the  system’s  brightness  outside  eclipsei  so 
the  intrinsic  variation  of  the  fainter  star  was  not  uniform  over  its  sur- 
face. A similar,  but  a much  smaller,  difference  existed  between  the 
other  observed  primary  minima.  The  observed  difference  in  depth  of 
primary  could  be  explained  by  a change  in  the  radius  of  the  fainter  star. 
Unfortunately,  there  were  insufficient  data  to  determine  if  this  was  the 
case.  Another  possible  explanation  would  have  been  a change  In  the  color 
of  the  fainter  star.  The  data  available  to  Pepper  indicated  a change  in 
cp1ar  of  the  pnoper  sense  but  of  insufficient  site  to  account  for  the 
observed  difference.  Also,  since  the  variations  outside  eclipse  had  Uie 
same  color  as  the  fainter  star  it  was  very  reasonable  to  expect  this  star 
to  be  the  source  of  Che  variation.  These  two  effects  (variation  in  Che 
depth  of  eclipse  and  the  variations  outside  eclipse)  could  have  been  due 
to  a combination  of  pulsation  and  "spottedness"  of  the  fainter  star, 
according  to  Popper.  He  pointed  out  that  the  depth  of  secondary  eclipse 
as  observed  by  Keller  and  Limber  (19S1)  was  greater  than  expected  from 
their  values  of  k and  the  depth  of  primary  minimum.  This  led  him  to 
believe  that  there  was  a large  uncertainty  in  the  value  of  k.  Popper 
also  found  a discrepancy  in  the  ratio  of  surface  brightness  as  calcu- 
lated from  the  spectral  types  and  that  obtained  from  the  light  curve.  He 
noted,  hewever,  that  the  discrepancy  was  probably  not  as  large  as  that 
found  by  Keller  and  Linier.  For  reference,  the  results  of  Popper  are 
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Clifsarl  and  Lacona  (1965)  oKerved  9S  CVn  from  tlie  S.  Agata  LI 
eattlatl  station  using  a 30  on  (13  Inch]  Cassegrain  telescope  In  1963 
and  1964.  Their  effective  wavelength  was  stated  to  be  5150  A‘,  which  is 
slightly  bluer  than  that  for  the  standard  V filter.  The  bandwidth  was 


not  given,  but  the  6HI  6366A  ohotomultiplier  tube  was  filtered  with  a 
yellow  Galileo  6 1-26  filter.  They  used  the  sane  star  as  Popper  for  a 
ccanparison  (M011477a  ■ BOwSS"  2420).  However,  they  did  not  correct  the 
observations  for  the  effects  of  aOnospherlc  extinction,  because  they 
considered  the  closeness  of  the  comparison  to  the  variable  sufficiently 
small  to  make  correction  unnecessary.  After  rectification  of  the  light 
curve  they  solved  the  1964  observations  to  obtain  the  following  elertents 
Initial  epoch  [Julian  Date]  > 2433467. 12B2, 


They  found  the  period  to  be  decreasing  in  the  sane  sense  as 
Keller  and  Limber  (1961)  had  found,  but  not  by  the  amount  predicted 
by  observations  obtained  prior  to  1936.  They  concluded  that  the  systen 
was  semi-detached  and  that  the  material  one  would  expect  in  such  a 
systen  was  leaving  the  larger  star.  Thus  a cloud  would  be  formed  which 
would  cause  the  ancnelous  light  variations  in  the  light  curve.  This 
does  not  agree  with  subseouent  research. 

Plavec  (1967)  studieo  the  properties  of  Algol-like  close 
binaries  with  double-line  spectra.  R5  CVn  was  among  the  stars  studied 
and  was  found  to  be  deteched  by  a significant  amount.  Several  other 


stars  were  noticed  to  liave  sintllar  properties  wMcn  were  different  from 
genuine  Algol  Oinaries.  RS  CVn,  UW  Dra,  Z Her,  AR  Lac.  and  SZ  Psc  were 
considered  by  Plavec  to  forsi  a well'defined  sub-group.  It  is  significant 
tnat  except  for  SZ  Psc  all  of  these  stars  were  in  Hiltner's  list  of  13 
systenB  showing  Ca  II  eotission. 

Catalano  and  Rodono  (IddZ)  presented  additional  photoelectric 
observations  of  RS  CVn  taken  with  the  same  equicment  Chisari  and  Lacona 
had  used  in  1963  and  1964.  However,  it  may  be  of  importance  to  note 
that  the  observations  were  not  all  taken  at  the  same  site.  The  data  for 
1965  were  apparently  obtained  from  the  same  place  as  those  of  Chisari 
and  Lacona,  bint  the  equipment  was  then  moved  to  Serra  La  Have.  The 
subsequent  data  taken  in  1966  were  from  this  site  using  the  same  tele- 
scope and  equipment.  All  of  the  observations  used  Che  same  comparison 
star  as  Popper,  Chisari  and  Lacona,  and  the  present  writer;  namely 
BW6°  Z420-  They  also  observed  a check  star  (BIH3S*  ZA22).  which 
showed,  contrary  to  Popper's  earlier  suspicions,  that  Me  comparison  star 
3Cif35^  2420  is  of  constant  light.  They  made  nn  attenpc  at  a solution,  in 
the  classical  sense,  to  the  new  observations,  but  instead  they  pointed 
out  an  important  fact  about  the  distortion  wave  in  RS  CVn. 

The  distortion  wave  was  found  to  shift  with  respect  to  the 
eclipses  in  the  direction  of  decreasing  phase,  and  completed  one  Qcle 
In  2400  orbital  periods.  To  explain  this  phenomenon  a model  was 
developed  in  which  a ring  of  material  circled  the  hotter  star.  This 
was  similar  to  the  proposal  of  Chisari  and  Lacona,  but  the  model 
developed  by  Catalano  and  Rodono  was  much  more  detailed.  In  this  mode! 
a non-uniform  ring  of  material  was  located  around  the  smaller  star  of 
the  system.  They  assumed  that  Che  material  would 


make  Xeplarian 


around  Oiia  star  and  -era  Uierefera  able  to  compute  the  distance 
beween  the  star  and  the  ring.  The  period  of  the  quasi-sinusoldal 
distortion  »a»e  represented  one  orbit  of  a particle  in  the  ring.  The 
results  required  that  the  material  be  outside  the  limiting  Boche  lobe 
for  the  star.  In  this  case  the  ring  »ould  be  unstable,  but  they  proposed 
that  the  material  lost  from  the  system  be  replaced  by  material  from  the 
larger  star.  If  the  primary  star  is  inclined  to  the  orbital  plane,  its 
axis  would  be  expected  to  precass  with  a period  which  they  computed  to 
be  about  the  same  as  that  of  the  retrograde  period  of  the  distortion 
wave.  The  quas1-sinuso1dal  shape  would  then  be  due  to  the  non-uhlform 
distribution  of  roterlal  In  the  ring,  but  there  was  no  exianatioh  for 
the  persistence  of  this  particular  distribution  over  the  time  scale  of 
the  observations  (1963  to  1966).  The  1965  light  curve  showed  no  appre- 
ciable variations  beyond  those  alreedy  known,  but  the  data  for  1966  In- 
dicated that  the  depth  of  primary  eclipse  was  not  the  same  for  all 
nights  during  which  this  portion  of  the  light  curve  was  observed.  The 
phenomenon  was  earlier  reported  by  Keller  and  Limber  and  by  Popper. 
Whether  thia  was  real  or  a result  of  the  new  observatory  site.  In  tne 
case  of  Catalano  and  Rodono's  data,  is  not  known. 

nelson  and  Duckworth  (1968)  retorted,  in  an  abstract,  on  the 
observations  they  had  obtained  of  RS  CVn  from  1966  to  1967.  Their 
observations  were  not  given,  but  they  reported  that  the  light  curves 
show  variations  in  the  depth  of  primary  and  secondary  eclipses  of  up 
to  a duarter  of  a magnitude  In  the  violet,  and  similar  but  smaller 
variations  In  the  other  colors.  They  also  reported  that  the  data  showed 
changing  shapes  for  the  eclipses  and  Hie  portions  of  the  light  curve 
outside  eclipse,  over  the  three  years  they  observed  the  system. 


Catalano  and  ftodono  (1966)  sunnarfzed  aaveral  aspects  of  their 
observations  of  RS  CVn  1n  a more  widely  accessfble  fonp  than  ttieir 
previous  publication.  First,  they  found  additional  evidence  for  the 
distortion  wave's  retrograde  shift.  From  the  observations  of  Keller 
and  Umber,  Popper,  Chisarl  and  lacona,  their  own  published  data  for 
1965  and  1966,  and  new  unpublished  data  for  1967  and  1963,  they  were 
able  to  more  accurately  determine  the  period  of  the  retrograde  shift. 
Unfortunately,  for  their  previous  model,  the  results  were  Incompatible 
with  the  theoretical  procession  rate  of  2100  orbital  periods.  They 
found  that  the  distortion  wave  retrograde  shift  period  was  of  Che  order 
of  800  orbital  periods  or  less.  The  new  data  for  1967  and  1968  were  ob- 


tained on  the  standard  UBV  system,  which  allowed  them  to  detennine  the 

and  outside  eclipse  varlat1o^ It  was  obvious  that  the  color  Index  In- 
creased with  decreasing  Intensity  of  the  distortion  wave.  Sinoe  the 
color  of  the  distortion  wave  became  bluer  as  It  became  fainter,  they 
suggested  that  the  Investigations  of  Mergentaler  (1950)  could  De  a 
passible  explanation.  They  reported  that  Mergentaler  proposed  that 
systems  which  snow  such  a color  change  can  be  explained  by  a gas  of 
negative  hydrogen  ions  of  different  optical  tnicknesses.  Catalano  and 
Rodono  noted  tmt  this  would  require  that  the  gas  be  In  an  eguilibriun 
configuration  In  order  to  maintain  the  shape  of  the  distortion  wave  as 
they  had  observed  It.  Variations  in  the  depth  of  primary  minimum  had 
been  observed  in  the  past,  but  their  more  consistent  and  more  extensive 
observations  were  able  to  confini  this  ohencmenonand  to  Show  that  it  was 
correlated  with  the  position  of  the  distortion  wave.  Since  the  primary 
ecHpse  was  total,  only  the  secondary  (larger,  cooler  star)  of  the  system 


proDably  responsible 


tbe  distortion  wave  phenomenon.  It  is  not  dear  to  the  oresent  writer 
if  Catalano  and  Rodono  were  referring  to  the  secular  C800  orbital  period) 
shift  of  the  distortion  wave  across  the  primary  eclipse  as  the  source  of 
variations  in  its  depth  or  to  some  sort  of  flaring  activity  on  the  larger 
star.  In  the  former  case  it  would  be  difficult  to  explain  changes  in 
depth  of  primary  eclipse  which  have  been  observed  on  much  shorter  tine 
scales.  Nevertheiess.  their  conciusion  was  in  complete  agreement  with 
the  eariier  conclusions  of  Payne-GaposchXin  {1939],  Kelier  and  Limber 
(1951),  and  Popper  (1961)  regarding  the  source  of  the  light  variations 
(in  eclipse  and/or  out  of  eclipse). 

The  period  variations  had  been  well  Known  previously,  but  there 
had  been  significant  evidence  that  Che  variations  were  due  to  apsids! 
motion  (Keller  and  Limber  1951;  and  Plavec  and  Smetanova  1959).  However, 
more  recent  studies  of  Plavec,  Smetanova  and  Pekny  (1961)  have  cast 
considerable  (fcubt  on  the  appropriateness  of  apsidal  nsition  to  explain 
the  period  variations.  Catalano  and  Rodono  showed  that  the  existence  of 

highly  questionable.  In  fact  they  stated,  ".  . . the  displacement  of 
secondary  minimum  resulting  from  our  observations  is  incompatibly  smaller 
than  the  orbital  period  variations"  (p,  440). 

Finally,  the  Catania  astronomers  concluded  the  article  with  some 
notes  on  the  spectroscopic  peculiarities  In  RS  CVn.  They  had  observed 
the  same  emission  lines  of  N^,  and  Ca  II  as  Hiltner  (1947),  Joy  (1930), 
and  Popper  (1961).  Unlike  KUtner  they  were  unable  to  notice  the  disap- 
pearance of  Che  H and  K lines  at 


secondary  eclipse. 


In  a tnaoretlcal  DresenUtlon,  HalT  11972)  discussed  an  alterna- 
tive fflodel  for  R2  CVn.  This  model  atteeteted  to  explain  the  following 
characteristics  of  Che  system: 

1.  The  existence  of  the  distortion  wave 

2.  The  migration  of  Che  wave  to  decreasing  phase 

3.  The  ancxnalous  depth  of  secondary  eclipse 

A.  The  non-un1fonn  migration  race  of  the  distortion  wave 

5.  The  variaOle  depth  of  primary  minimum  as  a function  of  time 

6.  The  displacement  of  secondary  minimum  from  a position  mid-way 
between  adjacent  pnlmary  minima 

7.  The  variable  amplitude  of  the  distortion  wave 
8-  The  changes  In  orbital  period 

The  model  that  was  proposed  Oy  Ha11 , Oy  analogy  with  the  sun,  was 
that  of  spots  on  the  cooler  component  of  the  system.  It  was  well  known 
that  the  sun  has  dark  spots  on  its  surface  from  time  to  time-  Starting 
from  this.  Hall  felt  that  It  was  reasonable  to  expect  that  other  stars 
would  have  similar  spots.  In  the  case  of  AS  CVn,  he  anticipated  tnat  Che 
spot  or  spots  would  cover  between  301  and  601  of  the  facing  stellar  disk 
of  the  XQ  star.  This  spotted  region  was,  for  whatever  reason,  confined 
to  one  hemlspnere  and  to  the  eguatarlai  region  between  about  r 30*^,  Hall 
reasoned  that  the  star  would  be  rotating  differentially,  again  by  analogy 
with  the  sun.  The  eguatorlal  region  would  rotate  faster  than  the  polar 
regions,  as  Is  the  case  In  the  sun  {which  also  has  a convective  envelope). 
The  stars  of  the  system  were  expected  to  be  in  synchronous  rotation,  but 
Che  latitude  at  which  the  rotation  was  synchronized  with  Che  or01t  was 
expected  to  be  of  the  order  of  * 30*^.  Again  relying  on  the  analogy  with 
the  sun,  he  assLUied  that  the  soot  activity  would  be  periodic,  with  a 


cycle  of  1800  orMtal  periods.  Fortner,  ne  assumed  that  tne  spots  on 
me  cooler  component  of  FIS  CVn  would  form  at  hipner  latitudes  at  the 
beginning  of  each  spot  cycle,  and  proceed  to  form  at  lower  and  lower  lati- 
tudes as  the  cycle  progressed. 

With  this  Timdel  the  existence  of  the  distortion  wave  is  easily 
explained.  Oifferential  rotation  of  the  star  explains  why  the  distortion 
wave  moves  to  decreasing  phase.  Since  Che  surface  of  the  cooler  star  is 
not  of  unifonn  brightness,  the  depth  of  primary  eclipse  will  vary  depend- 
ing on  Che  location  of  the  spotted  henisphere  at  the  time  of  this  eclipse. 
Furthermore,  the  location  of  the  spotted  region  can  cause  the  location  of 
the  miniimin  during  the  secondary  eclipse  Co  be  shifted  to  either  side  of 
phase  0.5.  Since  the  model  assianes  that  the  spot  activity  is  periodic  In 
a fashion  similar  to  the  sunspot  cycle,  and  that  the  star  rotates  differen- 
tially, then  the  migration  rate  of  the  distor-tion  wave  will  vary  depending 
on  the  latitude  of  the  spotted  region  and  the  amplitude  of  the  distortion 
wave  will  vary  in  a similar  fashion. 

Nall  presented  very  convincing  evidence  to  support  this  mode!. 

The  spotted  region  would  only  have  to  be  1000°i:  cooler  than  the  surround- 
ing photosphere  In  order  to  account  for  the  typical  amplitude  of  the 
distortion  wave.  This  was  of  the  proper  value  to  also  explain  the  anoma- 
lous depth  of  secondary  eclipse,  if  it  was  assumed  that  the  spotted  re- 
gion was  eclipsed  at  secondary  minimian.  Kail  noted  that  this  implies  that 
as  the  spotted  region  migrates,  due  to  differential  rotation,  around  the 
star,  the  depth  of  secondary  eclipse  should  vary  accordingly.  This  has 
not  been  confirmed. 

Catalano  and  Rodono  (1968)  bed  shown  the  correlation  between  the 
of  the  distortion  wave  and  the  depth  of  primary  minimum.  Hell 's 


poslti on 


(TOdel  provided  an  explanation  for  tMs.  The  deepest  primary  minleia 
occurred  durinp  times  when  differential  rotation  brought  the  spotted 
region  to  the  hemisphere  facing  the  earth.  The  dlsplacefnent  of  secon- 
dary eclipse  from  the  nld-polnt  of  the  light  curve  had  been  explained  by 
Catalano  and  Rodono  Cl^B)  as  a consequence  of  the  migration  of  the  dis- 
tortion wave.  Ball's  model  clarified  this  1n  terms  of  the  location  of 
the  spotted  region  during  the  secondary  eclipse.  If  Che  spotted  region 
was  located  so  that  1C  faced  the  earth  at  quadrature,  then  the  loss  of 
light  during  secondary  eclipse  would  be  asymmetric.  From  a plot  of  the 
available  estimates  of  the  amplitude  of  the  distortion  wave  as  a function 
of  epoch.  Hall  was  able  to  represent  the  variation  with  the  periodic 
function 

4V(max-m1n)  ■ -07l2  - (fo?  sin  ((E*4S0)/18Q0) . 
fie  used  the  linear  ephemeris  of  5chnel1er  (1926)  to  calculate  Che  value 
of  E.  Kota  that  this  function  can  be  used  to  predict  the  value  of  the 
amplitude  of  the  distortion  wave  at  any  time.  Rail  also  presented 
another  graph  which  allowed  the  prediction  of  the  distortion  wave  migra- 
tion rate.  In  Fig.  2 of  his  paper  he  plotted  the  orbital  phase  af  the 
minimum  of  the  distortion  wave  as  a function  of  epoch  (E).  From  this 
graph  one  was  able  to  see  Chat  the  discontinuities  In  the  migration  rate 
occur  at  the  minimum  of  tha  spot  cycle.  This  was  the  result  expected 
from  the  analogy  with  the  sun  and  differential  rotation.  Arnold  and  Hall 
(1973)  later  made  corrections  Co  the  model  without  changing  the  substance 
of  Che  results. 

In  the  remainder  of  Che  paper  He11  attempted  to  explain  the 
period  verletions  as  the  result  of  mess  ejections  from  the  spotted  active 
region.  The  required  mass  loss  rate  was  rather  high  (of  the  order  10"” 


$oUr  iD9$se&  per  year).  The  CQnduslon  he  reached  on  the  evolutionary 
status  of  RS  CVn  was  not  decisive,  but  he  considered  the  star  to  be  most 
likely  1n  a pre-nain-sequence  state  of  evolution. 

BS  CVn  Binaries 

The  early  stages  in  the  development  of  the  definition  of  this 
group  of  binary  stars  have  been  discussed  In  the  preceeeding  paragraphs. 
Certainly,  the  research  results  of  Struve  (1946)  and  KUtner  (1947)  were 
the  first  clues  to  the  existence  of  a group  of  stars  which  are  different 
from  the  classical  binary  and  significantly  different  from  other  peculiar 
systems.  The  research  effort  over  the  subsequent  30  years  (firatton  1950; 
Plavec  1967;  Popper  1967,  1970;  Oliver  1971,  1973,  1974;  and  Kail  I97B) 
slowly  developed  a clearer  picture  of  the  physical  and  observational 
characteristics  of  these  stars  which  separates  them  from  any  other  group 
which  has  been  classified. 


Popper  (1967,  1970}  listed  approximately  20  systems  which  he 
believed  to  be  of  special  interest  because  of  their  spectroscopic  charac* 
terlstics,  light  variations,  and  unique  placement  in  the  K-ft  diagram. 
Oliver  (1974)  was  the  first  to  undertake  a comprehensive  study  of  these 
stars  as  a group.  This  research  led  diver  (1971.  1973,  1974)  to  list 
the  criteria  for  inclusion  in  the  group  and  to  propose  a list  of  stars 
which  would  therefore  be  members.  To  date,  the  best  review  of  the  RS  CVn 
binaries  Is  the  recent  paper  by  Hall  (1976).  In  the  paragraphs  to  follow, 
the  current  understanding  of  the  RS  CVn  phenomenon  will  Pe  discussed. 

Since  the  most  characteristic  oOservationel  property  of  the  group  is 
strong  H and  K emission,  the  discussion  will  begin  with  the  spectroscopic 
characteristics.  Variations  in  the  observed  light  of  the  systems  are  an 
observational  property  of  many  of  the  systems  in  the  group.  The  discus- 


Sion  of  tliese  varlstions  will  be  followed  by  a short  report  on  ultraviolet 
and  infrared  e«ess.  Period  usriatiors  have  also  been  noticed  in  many  of 
the  RS  CVn  binaries;  the  short  and  long-term  variations  will  be  reviewed. 
The  recent  discovery  of  radio  and  soft  X-ray  emission  from  several  of  the 
RS  CVn  binaries  has  been  very  exciting.  The  report  here  will  attemot  to 
swnarize  these  observations  and  the  resulting  implications.  Finally, 
this  section  will  conclude  with  a review  of  the  proposed  evolutionary 
status  of  these  stars.  However,  before  beginning,  it  will  be  useful  to 
state  the  properties  of  RS  CVn  binaries  which  Hall  (1876)  proposed  as 
the  defining  characteristics.  They  are  as  follows; 

1.  Orbital  period  between  1 day  and  2 weeks 

2.  Strong  K and  K emission  seen  outside  eclipse 

The  list  of  manhers  in  this  group  which  has  been  developed  over  the  last 
few  years  is  given  in  Table  3.  This  list  is  derived  from  Hall's  (1976) 
review  caper  and  new  discoveries. 

Spectroscopic  characteristics 

The  existence  of  strong  K and  K emission  In  the  RS  CVn  binaries 
has  been  well  docLotenCed  by  Struve  (1946),  Hlltner  (1947),  Eggen  (1965), 
Popper  (1967,  1970),  Oliver  (1974)  and  Ueiler  (1978a).  The  emission  is 
strong  in  the  sense  that  it  is  usually  significantly  above  the  continuum. 
In  no  case  is  this  H and  X emission  to  be  confused  with  the  normal  emis- 
sion core  reversal  in  late-type  single  stars  which  gives  rise  to  the 
H1l50n-8appu  effect. 

Struve  (1946)  studied  the  variatioh  in  the  strength  of  the  K and 
K lines  as  a function  of  phase  angle  in  several  stars.  It  was  his  con- 
clusion that  the  maximum  strength  was  at  the  quadrStune  points  (0?25  and 
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0?75).  Oliver  (1974)  suimtarleed  the  existing  data  by  saying  that  the 
enrission  was  usually  at  a ninlmun  during  the  secondary  ecTiDse,  but  not 
in  all  cases,  nor  at  a11  times.  Hlltner  (1947)  had  found  a disappearance 
of  the  N and  K emission  lines  at  secondary  eclipse  in  H5  CVn.  This  led 
him  to  Oelleve  that  the  snission  originated  from  localized  sources  at 
the  BXtrenities  of  a tidal  bulge  along  the  line  Joining  the  two  stars. 
Struve  (1946)  had  earlier  offered  a similar  explanation  for  the  origin 
of  the  emission  in  similar  cases.  Catalano  and  Rodono  (1966)  were 
unable  to  confirm  the  disappearance  of  the  Ca  II  onlssion  lines  at 
secondary  mininnmi,  thus  raising  some  doubt  as  to  the  validity  of 
Kiltner's  tidal  bulge  hypothesis.  Here  recent  studies  by  Weiler  (1975a, 
197Ba)find  that  the  emission  is  associated  with  the  lerger  star,  but  the 
emission  Is  not  localized.  Furthermore,  Vfeiler  does  not  find  the 
strengthening  at  the  guadriture  points  or  the  tendency  for  the  emission 
to  weaken  at  secondary  minimum. 

The  existence  of  Ca  II  emission  Is  Indicative  of  chromospheric 
activity  (Toung  end  Konlges  1977;  Weller  1975a;  NaftHan  and  Drake  1977; 
and  Weller  1976a).  This  eii1sslon  must  originate  frcn  the  larger  (cooler) 
ccknponent  because  the  eoilssian  lines  have  velocities  which  agree  well 
with  the  absorption  lines  for  this  component  (Popper  1951;  Oliver  1974). 
In  addition,  the  width  of  the  emission  lines  Is  of  the  proper  amount  to 
be  explaioed  by  the  synchronous  rotation  of  the  larger  star  (Popper 


The  emission  lines  of  Ca  II  are  not  the  only  emission  lines  which 
nave  been  reported.  Weller  (1978a)stud1ed  6 RS  CVn  systems  and  found 
to  be  in  mission  In  all  6.  Naftllan  (1975)  found  to  be  in  enis- 


(weakly) 


associated  with  the  cooler  star  by  Kaftilan  and  Drahe  (1977).  The  H 
eoilssion  In  the  six  systens  studied  by  Weiler  (197Sa,  1978a)  was  variable 
to  a significant  degree  In  four  of  them.  Similar  emission  in  the  H Tine 
has  been  reported  by  8app  (1978)  for  Che  systems  H9  1899  and  UX  Ari. 
Hcxever,  it  was  noted  that  these  were  the  only  two  stars  from  a survey 
of  about  30  stars  which  shoved  persistent  H emission.  Three  other  stars 
showed  sporadic  H emission  at  the  resolution  used.  Weller  (1978b) 
points  out  that  at  higher  resolution  It  Is  possible  that  all  RS  CVn 
binaries  may  show  some  filling  in  of  the  K Tine  by  emission. 

The  H and  K and  X emission  in  several  RS  CVn  binaries  has  shown 
correlation  with  the  tfiese  position  of  the  distortion  wave  and  the  phase 
of  greatest  emission  (8opp  1978;  Herbst  1973;  and  Weiler  197S6,  1979a). 

In  additien  Weiler  finds  the  distortion  wave  more  pronounced  the  greater 
Che  correlation.  This  1s  strongly  suggestive  of  a physical  relationship 
between  the  causes  of  Che  two  phenomene.  Weller  (ig78a)  suggests  that 
the  emission  is  due  Co  an  active  photosphere  idrich  in  turn  drives  an 
active  chromps^ere  in  the  region  of  the  photospheric  activity.  Young 
and  Xoniges  (1977)  attribute  this  to  tidal  caupling  which  increases  the 
density  scale  height  and  thereby  changes  the  acaustic.wave  power  spectnim. 
These  Investigators  feel  that  the  resulta  support  Hal! 's  ( 197Z)  model 
for  a spotted  surface  on  the  cooler  component  in  RS  CVn.  and  possibly 
other  stars  with  similar  characteristics. 

Oliver  (1974)  found  a variation  in  the  absorption  line  strengths 
but  the  correlation  with  Tight  curve  variations  was  weak. 

The  abundance  of  metels  In  binary  systems  was  studied  by  Miner 
(1966).  In  the  systems  he  studied,  an  underabundance  of  metals  was 
found  in  eoUpsing  binaries.  The  urtderabundance  does  not  appear  to  be 
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as  great  in  the  3 RS  CVn  systems  which  were  inclLided  in  the  project  as 
far  the  ether  eclipsing  systems.  The  seb-giant  components  of  12  Algol* 
like  binaries  were  observed  by  Hall  (1967)  and  foirnd  to  have  an  under- 
abundance  Of  heavy  elenents.  NaftHan  (1975)  found  the  cooler  stars  in 
RS  CVn  and  RH  UMa  to  be  underabundant  in  iron  and  chromium,  which  he 
noted  is  in  contradiction  to  a pre-main-sequence  state  of  evolution. 

The  cooler  component  of  Aft  Lac  was  also  found  to  be  deficient  in  metals 
by  Naftilan  and  Drake  (1977),  Anderson  and  Popoer  (1975)  concluded  that 
from  the  observed  radius  and  temperature  of  the  stars  in  TV  Pyx  there  is 
an  indication  of  an  abundance  (or  opacity)  anoinaly  in  both  stars.  This 
would  seen  to  indicate  that  RS  CVn  binaries  are  not  very  young.  Haftilan 
and  Drake  (1977)  have  concluded  from  their  research  that  Che  anomalously 
low  metal  abundance  is  not  due  to  free-free  emission  filling-in  the  metal 
lines  in  AR  Lac,  nor  is  it  due  to  a circumstellar  shell.  Recent  soft 
X-ray  emission  from  UX  AH  has  led  Halter,  Charles,  and  Bowyer  (1978)  to 
conclude  that  the  heavy  element  depletion  in  this  system  is  real,  because 
they  do  not  observe  an  Fe  emission  line  at  0.8S  KeV  as  was  seen  in  the 
X-ray  spectrun  of  Capella, 

Lithium  (Li)  has  not  been  detected  in  several  RS  CVn  binaries; 
this  supports  the  view  that  these  stars  are  not  very  young  (Naftilan  and 
Drake  1977;  Young  and  Koniges  1977;  and  Conti  1967).  Dr,  the  other  hand, 
Rucinski  (1977)  has  found  a weak  Li  line  in  HD  224086.  This  system  was 
recently  added  by  Hall  (1978)  to  the  list  of  RS  CVn  type  binaries. 

As  has  been  mentioned  previously,  a correlation  between  chrono- 
spneric  Ca  11  emission  and  tidal  coupling  has  been  found  by  Sratton 
(19SD),  Young  and  Xoniges  (1977),  and  GUbocki  and  Stawikowskf  (1977). 
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)oimg  and  Koniges  believe  that  tnis  is  evidence  tnat  when  a star's  outer 
atmosphere  is  affected  by  gravitational  attraction  of  a close  companion 
Che  result  is  a change  in  the  density  scale  height.  This  shifts  the 
acoustic-wave  power  spectrum  to  lower  frequencies,  which  results  in  a 
large  increase  in  the  dissipation  of  energy  into  the  chromosphere. 

Gratton  Delieved  that  the  tidal  distortion  would  produce  localised 
regions  of  emission  in  accordance  with  Hiltner's  model.  However,  the 
latter  interpretation  was  eicluded  by  Struve  (1948)  lAen  he  was  unable 
to  find  the  expected  degree  of  linear  polarization  from  such  a localized 
region.  Weiler  (1978a)was  also  unable  to  find  the  variation  in  the  H 
ahd  K emission  lines  as  a function  of  phase  that  would  be  expected  from 
the  tidally  distorted  bulge.  He  concluded  that  the  emission  may  well  be 
somewhat  localized.  Out  not  at  the  extremities  of  a tidal  bulge  formed 
along  the  line  Joining  the  two  stars.  Glebccki  and  Stawitowski  believe 
that  their  study  of  H and  k emission  supports  the  tidal  bulge  hypothesis, 
while  Heiler  feels  that  his  research  clearly  points  to  an  active  chromo- 
sphere, Time  may  resolve  the  dispute. 

Photometric  anomalies 

The  light  curves  of  many  (if  not  ell)  of  the  RS  CVn  binaries  show 
various  anomalies.  The  most  conmon  and  striking  light  curve  anomaly  is 
a quasi'sinusoidal  wave  that  distorts  the  normal  eclipse  phenomenon  as 
can  be  seen  in  Figure  1,  This  wave-like  distortion  has  been  referred  to 
as  the  "distortion  wave"  by  Oliver  (1974);  this  term  will  be  used 
throughout  this  dissertation.  There  are,  however,  other  phenomena  which 
are  abnormal  in  the  light  curves  of  these  binaries.  Several  of  the  sys- 
tems show  season-to-season  changes  in  the  system  brightness,  while  a few 
have  exhibited  increases  in  brightness  over  very  short  time  scales.  The 
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short  tine  scale  variations  have  been  referred  to  as  flare-1ike  activity 
by  some  observers.  These  Irregular  changes  in  the  light  of  the  system 
cannot  be  explained  by  changing  parameters  of  the  distortion  wave,  and 
are  therefore  discussed  as  separate  (but  perhaps  related)  phencmena. 
aesldes  the  variations  In  the  light  of  the  systan.  the  light  curves  of 
these  binaries  have  two  additional  anetnalies  which  warrant  inentlon. 
first,  the  depth  of  the  primary  eclipse  in  several  systems  has  shown 
vanlatlpns,  which  are  not  associated  with  changes  in  system  brightness. 
Second,  the  secondary  eclipse  for  many  systans  has  shown  variations  In 
Its  position  anO  depth. 

The  distortion  wave  in  fis  CVn  was  first  noted  by  Sltlerly  (1930) 
as  a difference  In  the  heights  of  the  maxinia.  This  observation  was  noted 
by  several  researchers  before  the  true  picture  was  developed  by  the 
Itelian  astronomers  at  Catania.  Following  the  advice  of  Popper  (1961) 
the  Catania  astronomers  (Chlsarl,  Lacona,  Catalano,  and  RoOono)  started 
a long-term  program  of  photoelectric  observations  of  fiS  CVn,  TMs  work 
produced  the  remarkable  discovery  that  the  difference  In  the  heights  of 
the  maxima  was  changing,  and  the  changes  could  be  explained  by  a mlgnating 
quasi-s1nuso10il  wave  superimposed  on  the  light  curve.  Chlsarl  and 
Lacona  (1965)  and  Catalano  and  Rodono  (1967,  1968)  showed  very  convincing 
evidence  for  the  existence  of  this  wave,  and  for  the  slow  migration  to 
decreasing  phase  angle  from  season-to-season.  They  referred  to  the  slow 
shift  of  the  distortion  wave  to  decreasing  phase  as  a “retrograde  nigre- 
tion.'’  Oliver  (1974)  has  shown  that  this  property  is  very  common  in  the 
RS  CVn  binaries  which  he  studied.  It  Is  ccmnonly  held  that  the  existence 
of  the  distortion  wave  and  the  retrograde  m1gret1on  are  properties  that 
ell  RS  CVn  binaries  shere.  The  evidence  Is  not  complete  for  ell  the 


systems  cansidsred  to  be  of  ttie  RS  CVn  oinery  type,  but  many  systems  tieve 
shown  te1s  phenomenon  {8opp,  Espenak,  Ka11,  LandU,  love11,  and  Reucroft 
1977;  Oliver  1974;  Hall  1972,  1977;  Popper  1974;  QiantUss  1976;  Hilone 
1969;  iielspn  and  Duckworth  1969;  Nell,  Henry,  Ourke,  and  Mullins  1977; 
Hall,  Hpntle,  and  Atkins  1975,  Blanco  and  Catalano  1970;  Landis,  Lovell, 
Hell,  Henry,  and  Renner  1978;  Neernshaw  end  Oliver  1977;  and  Rucinski 
1977).  Ikany  of  the  above  Investigators  reported  migration  rate  changes, 
and  variations  in  the  amplitude  of  the  distortion  wave  from  season-co- 
season.  Perhaps  significantly,  Eaton  (1977)  has  found  the  distortion  wave 
to  be  moving  to  increasing  phase  in  the  system  52  Psc.  The  retrograde 
migration  in  SS  8oo  has  apparently  reversed  itself  according  to  Hall 
(1978).  The  color  dependence  of  the  distortion  wave  amplitude  is  very 
interesting.  For  nost  of  the  systems  (reported  in  the  aforementioned 
publications)  the  amplitude  is  greater  in  the  V band  than  in  the  B band. 
However,  there  are  notable  ekceptions  to  this:  M Ore  (Oliver  19741, 

AT  Lac  (Hilone  1968),  UV  Psc  (SaOik  1978),  and  SZ  Psc  (Eston  1977). 

In  sdOition  to  variations  in  the  amplitude  of  the  distortion 
wave  within  a system,  the  siae  of  the  amplitude  varies  a great  deal  from 
system  to  system.  The  largest  amplitude  appears  to  be  ahbut  0?2  for 
RS  CVn  (in  tne  V band),  and  the  smallest  appears  to  be  abbut  0?02  for 
several  other  systems.  Recently,  infornetlon  ss  to  the  actual  shape  of 
the  distortion  wive  has  become  available,  ftjcinski  (1977)  noted  an  Inter- 
esting change  in  the  shape  of  the  light  variations  for  HD  224085,  and  the 
present  writer  (Ludington  1976]  has  noted  changes  in  Che  detailed  shape  of 
the  distortion  wave  in  RS  CVn,  Hall,  Hontle,  and  Atkins  (1975)  noticed  a 
change  in  the  shape  and/or  amplitude  of  UX  Ari  from  their  UBV  and  JHKL 
photometry  in  1972.  The  significance  of  these  variations  is  not  clear, 
but  they  would  seem  to  point  to  dynamic  structures  on  the 
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Copper  U967]  anil  Oliver  <197<)  nave  noted  that  the  RS  CVn  sys- 
tems show  small  Irregular  variations  in  the  system  brightness.  Wod  (1946) 
was  perhaps  the  first  to  observe  these  irregular  variations  in  the  system 
AR  Lac.  Nora  recently,  Bopp  et  al.  (1977)  have  found  evidence  for  irreg- 
ular variations  of  iO'?B2  in  the  output  of  HR  1099.  Larger  variations 
have  been  discovered  in  the  brightness  of  UX  Ari  by  Hall  (1977)  and 
Landis  et  al.  (1978).  The  latter  group  also  detected  a flare-lihe  event 
of  JO  2443054.76  vdiich  lasted  for  approximately  one  hour,  and  Increased 
the  system  brightness  by  0?15. 

The  Irregular  variations  discussed  in  the  last  paragraph  refer 
to  the  total  system  brightness.  Light  variations  have  been  reported 
during  the  primary  eclipse  by  Hood  (1946)  and  Keller  and  Lindjer  (1951) 
for  AR  Lac  and  R5  CVn,  respectively.  Similar  changes  during  primary 
eclipse  have  been  noted  by  Popper  (1961)  for  RS  CVn:  Hilone  (1977)  for 
RT  Lac:  Nelson  and  Duchworth  (1968)  for  RS  CVn;  and  Oliver  (1974)  for 
SS  Boo,  HV  One,  RS  CVn,  I Her.  AR  Lac,  RT  Lac,  and  S2  Psc.  It  should 
be  noted  that  several  of  these  systems  have  total  (occultation)  primry 
eclioses.  Therefore,  if  the  variation  is  seen  during  this  total  phase, 
it  is  expected  that  the  variation  is  due  to  intrinsic  variability  of 
the  larger  (cooler)  aar  only.  In  the  oase  of  AR  Lac.  this  is  not  the 
case.  The  total  portion  of  primary  eclipse  was  seen  to  have  small  or 
no  variation  (Kron  1947).  Kron  found  that  the  variation  occurred duri ng 
the  ingresa  and  egress  branches  of  the  eclipse. 

Since  the  primary  eclipse  in  RS  CVn  is  total,  then  the  secondary 
eclipse  will  be  annular.  Thus,  since  the  larger  star  is  still  visible. 

It  is  not  surprising  that  this  eclipse  also  shows  anomalous  behavior. 


Jellap  and  Limber  0951)  found  that  frcm  ttia  depth  of  the  ecliose  the 
ratio  of  Intensities  was  much  less  than  would  be  expected  for  stars  of 
the  spectral  types  which  make  up  RS  CVn.  Popper  [1961]  found  the  sane 
Inconsistency  with  the  system,  but  also  found  the  discrepancy  to  be  less 
severe.  In  addition,  Keller  and  Linber  reported  a small  eccentricity 
from  the  displacement  of  secondary  eclipse,  neither  the  spectroscopic 
orbit  of  Joy  (1930),  nor  the  orbit  of  Popper  0961)  Showed  anything  but 
a circular  orbit.  The  distortion  wave  discovered  by  the  Catania  ob- 
servers suooHes  the  explanation  for  both  the  changing  depth  of  primary 
minimum  and  for  the  displacement  of  the  secondary  eclipse  from  the  0?S 
position.  The  variation  In  brightness  at  the  total  phase  of  primary 
eclipse  and  the  dlsoUcement  of  secondary  eclipse  are  correlated  with  the 
migration  of  the  distortion  wave.  Hall’s  OS72)  model  took  the  distortion 
wave  explanation  one  step  farther.  If  the  sootted  region  that  Kali  pre- 
dicts does  indeed  exist,  then  the  encmalous  depth  of  secondary  eclipse 
is  explained  by  the  variation  of  the  surface  brightness  bf  the  larger 
star.  Kail's  spot  model  also  explains  the  displacement  of  secondary 
eclipse  in  more  concrete  terms. 

Optical  and  radio  oolarltatlon 

The  polanxation  of  the  electromagnetic  radiation  from  SS  CVn 
binaries  was  first  searched  for  by  Struve  (IMS).  In  this  study  he 
attempted  to  locate  the  linear  polarization  which  would  result  from 
the  tidal  bulge  proposed  by  Hlltrer  (1947)  as  the  origin  of  the  Ca  II 
en1s5lon.  The  H and  K lines  of  dll  Lac  showed  no  evidence  of  linear 
polarization.  Struve  placed  an  upper  limit  of  1OT  on  the  polarization 
from  this  star.  In  a survey  of  unevolved  and  evolved  main-sequence 
biniry  systwis  Pfeiffer  and  Koch  (1977a)  found  that  the  unevoived 


systems  showed  no  significant  linear  polarization. 


evolved  binaries 


showed  intrinsic  linear  polarization  if  the  log  S was  greater  than  or 
approximately  eguai  to  1.0.  Here  the  *S’  represents  the  separation  of 
the  stellar  photospheres  in  units  of  solar  radii.  They  further  stated 
that  the  incidence  of  polarization  in  this  group  does  not  depend  on  nass- 
ratiOi  total  system  masSi  fractional  component  size,  rotational  velocity, 
the  stage  of  evolution  (Case  A on  Case  B},  or  the  detached  or  semi- 
detached state  of  the  system.  It  was  their  belief  that  this  was  evidence 
for  the  necessity  of  sufficient  "free"  volume  to  allow  more  material  and/ 
or  greater  asynioetry  to  produce  the  linear  polarization. 

Three  AS  CVn  binaries  have  been  reported  to  have  some  degree  of 
polarization.  RS  CVn  was  reported  to  display  variable  linear  polariza- 
tion by  Pfeiffer  and  toch  (1973).  The  variation  was  cyclic  with  orbital 
period  and  showed  random  short-term  variation.  This  was  suggestive  to 
them  of  electron  and/or  Rayleigh  scattering  in  an  active  circumstellar 
envelope.  Only  an  upper  limit  could  be  set  by  Pfeiffer  and  Koch  (1977b) 
for  the  linear  polarization  in  the  visual  band  for  HR  1099  (less  than 
0.02X).  At  radio  wavelengths  no  significant  linear  polarization  has  been 
detected  for  HR  1099  or  UK  Ari,  but  a significant  degree  of  circular 
polarization  has  been  detected  (Owen  et  al.  1976i  and  Spangler  1977). 

The  degree  of  circular  polarization  was  found  by  them  to  be  a function 
of  wavelength,  and  was  only  detected  during  outburst  activity.  This  was 
suggestive  of  microwave  bursts  like  those  observed  from  the  sun. 

Pfeiffer  (1978)  recently  reported  on  some  very  detailed  polariza- 
tion studies  of  RS  CVn.  He  concludes  that  the  polarization  is  due  to 
scattering  in  a large  doud  which  engulfs  the  system.  He  also  conments 
that  fiS  CVn  itself  is  the  only  RS  CVn  binary  in  which  he  has  been  able 
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to  detect  linearly  polarized  light,  as  of  the  time  of  his  presentation. 
Ultraviolet  and  Infrared  eaeess 

Oliver  (19741  has  done  the  most  contprehensive  study  of  ultra- 
violet (UV)  excess  in  the  RS  CVn  binaries.  It  was  his  conclusion  that 
an  ultraviolet  excess  uas  a coimon  oceurrencein  the  U-8  color  index  of 
the  cooler  conponent.  Since  fee  of  the  stars  have  UBV  data  available, 
it  is  conceivable  that  virtually  all  the  systems  have  a cooler  component 
ivith  a UV  excess.  Recently,  Rhombs  and  Fix  {1977}  presented  data  which 
confimed  very  clearly  the  UV  excess  in  RS  CVn,  AR  Lac,  and  UX  Ari.  They 
also  were  able  to  show  that  the  cooler  component  was  responsible  in  each 
case.  Hall  (1972J  suggested  that  the  excess  could  possibly  be  analogous 
to  the  excess  in  T Tauri  stars.  Rhombs  and  Fix,  on  the  other  hand,  con- 
cluded that  free-free  emission  from  a hot  circumstellar  gas  was  the  best 
explanation. 

Infrared  (Ifi]  excess,  liXe  the  ultraviolet  excess,  appears  to  be 
a coimion  cbaracten'stic  of  nearly  all  RS  CVn  binaries.  The  problem  of 
the  infrared  excess  has  been  discussed  by  several  authors  (Atkins  and 
Hall  1972;  Hall,  Hontle,  and  Atkins  1975;  Hall  1976;  and  Milone  1976a. 
1976b).  There  is  general  agreement  thit  at  least  the  majority  of  RS  CVn 
binaries  display  sme  degree  of  infrared  excess,  but  the  source  of  the 
excess  is  controversial.  Nilone  (1976b)  found  an  IR  excess  in  10  of  14 
systems  showing  light  curve  asyemetri es . which  were  chosen  at  random  for 
study.  Of  this  group.  6 systems  ere  classed  by  Hall  (1976)  as  RS  CVn 
binaries.  Three  bf  these  were  definitely  found  by  Milone  to  have  an 
infrared  excess.  He  cautioned,  however,  that  the  other  two  stars  (SS  Boo 
and  Mi  Ora)  may  well  snow  the  excess  because  the  spectral  type  and  light 
ratios  were  taken  from  the  literature  and  could  cause  a mis-evaluation  of 


the  e«ess.  Atkins  and  Hall  (1672)  found  an  infrared  excess  In  5 or  6 
systems  tor  which  they  had  sufficient  data  to  make  an  adequate  evaluation, 
Kiione  C1976a.bl  considers  the  18  excess  to  be  due  to  a circuaistellar  jas 

did  not  consider  a clroumstellar  cloud  a possible  explanation  for  the 
JHKL  pnotcmetry  they  had  obtained.  This  question  remains  unsettled,  but 
in  light  of  Che  polarization  data  of  Pfeiffer  {19781  it  would  seem  that 
a circuTstellar  cloud  may  be  the  correct  explanation  at  least  in  the 
single  case  of  RS  CVn  itself. 

Period  variations 

The  period  variations  In  many  of  the  systems  listed  by  Hall  (1976) 
are  very  pronounced.  He  noted  that  a linear  epheneris  in  some  cases  can 


quarter  of  an  orbital  cycle  in  only  10  years  or 


so.  Hall  suspects  that  virtually  all  of  the  systems  have  these  large 
oeriod  variations,  because  about  a third  are  known  to,  and  the  rest  have 
insufficient  data. 

Since  the  RS  CVn  binaries  are  detached  systems,  the  period  changes 
are  not  due  to  the  sane  mechanism  as  in  the  semi-detached  Algol-type 
binaries.  In  the  orevlous  section  on  the  history  of  RS  CVn  it  became 
apparent  that  the  period  changes  are  also  not  due  to  aosidal  motion  or 
light-time  effects  of  orbital  motion.  Therefore,  it  is  very  clear  that 
the  mechanism  responsible  is  rather  unusual. 

Flare-type  mass  ejection  from  one  hemisphere  hes  been  proposed 
as  an  explanation  for  the  short-term  perloO  changes  (Hall  1972;  Arnold 
and  Hall  1973s  Hall  1975b;  and  Hall  1976).  In  this  model  the  brighter 
hemisphere  Is  ejecting  matter  by  a High-velocity  impulse-type  mechanism. 
TMs  results  In  a correlation  between  the  observable  short-term  period 
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c)ian;es  and  the  observabU  pbasa  of  the  distortion  »ave  lalnlnum.  The 
mass  loss  rate  renoired  to  account  for  the  period  variations  is  rather 
large:  lO'Sl^year.  In  the  latest  version  of  this  model.  Hall  (1976) 
suggests  chat  the  effective  isonent  arm  is  some  Alfven  radius  so  that  the 
required  mass  loss  rate  is  reduced. 

Several  aspects  of  this  model  have  been  criticized,  Catalano  and 
Rodono  (1974)  objected  to  ejection  of  material  from  the  brighter  hemi- 
sphere on  the  grounds  that  one  would  expect  flare-type  mass  ejection  from 
the  active  spotted  hemisphere.  They  also  pointed  out  that  the  measurement 
of  the  tine  of  minimum  was  affected  by  the  deformation  of  primary  eclipse 
by  the  distortion  wave.  This  effect  would  have  the  same  sense  as  the 
correlation  of  distortion  wave  oriniiman  with  period  change.  Hall  (1975b) 

sense,  the  magnitude  of  the  0-C  variations  by  deformation  of  primary 
eclipse  were  25  times  toe  small. 

In  tne  case  of  AR  Lac  It  has  been  possible  to  use  s modification, 
for  this  specific  system,  of  the  model  to  explain  the  long-term  period 
changes  (Hall,  Richardson  and  Chaniliss  1976),  The  long-term  period 
variations  in  the  other  systems  have  not  been  adequately  explained. 

Further  research,  both  observational  and  theoretical,  is  greatly  to  be 

The  possibility  of  an  Enhanced  stellar  wind  in  RS  CYn  binaries 
has  been  proposed  by  Oliver  (1974),  Ulrich  end  Popper  (1974)  and  Popper 
and  Ulrich  (1977).  This  is  a conceivable  source  of  period  changes  and 
deserves  further  Investigation  as  pointed  out  by  Hall  (1976). 


EinIsslOT  at  radio  mveltnctht 


The  detection  of  an  as  CVn  binary  at  radio  wavelengths  by  Gibson 
and  Hjellming  (1974)  was  believed  by  then  at  first  to  indicate  mechanisms 
similar  to  those  for  known  Algol-11fce  flares.  The  spectral  index  of  the 
flare  was  definitely  non-themal  as  was  the  case  for  tne  Algol-like 
flares.  Later  wort  by  Gibson,  Hjellnring,  and  Owen  (1975)  showed  that 
the  flares  were  non-themal  in  IX  Ari  but  that  an  evolving  synchrotron 
source  was  a more  likely  mechanism  than  the  infall  of  natter. 

Gibson  et  a1.  (1975)  cofiipared  the  flare  frcm  UX  Ari  with  those 
of  UV  Ceti  stars,  and  found  that  there  was  little  similarity.  The  spec- 
tral index  for  a UX  Ari  flare  was  about  *0.2.  while  typical  values  for 
UV  Ceti  stars  range  from  -3  to  -5.  Even  at  the  end  of  the  flare  observed 
on  11.  12  August  1974  the  spectral  Index  of  UX  Ari  only  reached  -Q.6. 
Therefore,  tne  type  Of  flare  is  different  front  those  In  UV  Ceti  stars. 

This  was  supported  by  the  decay  time-scales  in  the  2 types  of 
flares.  Gibson  et  al.  (1975)  found  that  a UX  Ar1  flare  decayed  on  a 
time-scale  of  days.  This  was  about  100  to  200  times  the  decay  rate  of 
a flare  in  a typical  UV  Ceti  star. 

Circularly  polarized  emission  has  been  detected  at  radio  wave- 
lengths by  Spangler  (1977),  and  Owen,  Jones,  and  Gibson  (1975)  from 
2 AS  CVn  binaries  (UX  Ari  and  HR  1099).  The  level  of  oolariaation  was 
about  51  (for  UX  Ari)  to  20t  (for  HR  1099)  at  1400  MHz,  about  at  (for 
HR  1099)  at  S06S  MHz,  and  less  than  21  (for  HR  1099)  at  2G96  MHz.  The 
detection  of  circular  polarization  in  these  stars  contrasts  sharoly  with 
the  lack  of  circularly  polarized  emission  from  Algol. 

Ho  linear  polartzation  has  been  detected  et  radio  wavelengths 
from  any  RS  CVn  binary  at  the  time  of  this 


however  have  out  an  uooer  limit  of  2-3t  at  2695-8085  MHj  for  the  linearly 
DOlarited  enissioh  from  HD  1099. 

No  correlation  has  been  found  Secween  ootical  and  radio  variability 
as  reoarted  by  Chambliss  (1976)  and  Spangler  (1977).  This  is  not  a sur- 
prising result  if  the  flare-like  activity  is  similar  to  the  typical  solar 
microwave  flare.  In  these  solar  flares  the  energy  output  in  the  radio 
region  is  much  greater  than  in  the  optical  region. 

Radio  binaries  tend  to  be  of  late  spectral  type  (G  and  K),  and 
they  tend  to  be  above  the  main-sequence  according  to  Owen  and  Soangler 
(1977).  This  may  be  part  of  the  eaplanation  for  the  high  freoueney  of 
radio  emitters  among  RS  CVn  binaries.  On  the  other  hand,  and  perhaps 
more  likely,  the  statistics  may  Pe  biased  by  the  selection  of  RS  CVn 
binaries  as  candidates  for  radio  emission  surveys.  Wore  dataare  required 
to  establish  this  point  more  clearly  either  way. 

Owen  and  Spangler  (1977)  were  unable  to  detect  iny  change  in  the 
observed  flux  from  AR  Lac  at  4585  WHt  during  any  eclipse.  This  Is  a good 
indication,  as  they  pointed  out,  that  the  source  is  not  a compact  region 
between  the  Wo  stars.  They  computed  that  the  radius  and  brightness 
temperature  of  a source  located  in  this  region  would  be  limited  by 
r » J X lo”  on,  and  Tj^  . 4 x 10®  'K. 

Gyrosynchrotron  radiation  has  been  suggested  as  the  source  of  the 
radiation  [Gibson  et  al.  1975;  Owen  et  al,  1976;  Spangler  1977;  and 
Owen  and  Spangler  1977),  These  authors  have  found  that  this  mechanism 
explains  the  radio  flare  observations  better  than  other  alternatives  they 
have  considered.  It  accounts  for  the  observed  non-thermal  spectral  index, 
the  circularly  polarized  emission  during  flare  events,  and  the  absence  of 
linear  polarization  (If  reasonable  assumptions  are  made  about  the  physical 
properties  of  the  emitting  region). 


A magnetic  field  of  adoot  30  fiauss  and  an  electron  denalCy  of 
about  Z.i  a 10^^  cm  ^ would  be  indicated  by  Spangler'a  (1977)  calcula- 
tions. Spangler  et  al . (1977)  felt  that  this  imdel  core>Hnents  Mullen's 
(1974.  1976)  theory,  and  therefore  would  also  seen  to  support  the  star- 
spot  hypothesis  of  Hal!  (1972,  197$). 

Soft  7-ray  emission 

A oDst  exciting  discovery  has  been  the  recent  detection  of  soft 
t-ray  emission  (0.2  to  2.B  KeV)  from  mry  R5  CVn  binaries.  As  of  June 
1978,  Liller  (1978)  reported  that  II  of  the  RS  CVn  binaries  had  been 
identified  as  x-ray  sources.  This  compares  with  two  or  three  Identifica- 
tions reported  by  Charles  (1978)  in  April  1978.  This  discovery,  along 
Wia  the  previously  known  radio  emission,  gave  additional  credence  to  the 
solar  analogy  of  highly  aotive  regions  In  the  atnosohere  of  the  sub-giant 
component  as  pnoposed  by  hall  (1972,  1976). 

Cash,  8owyer,  Charles,  Lampton,  Sarmire,  and  Aiegler  (1978)  were 
able  to  fit  a solar  abundance  plasma  at  about  10^  to  the  soft  X-ray 
spectrum  of  Capellt.  This  star  was  Included  as  a related  system  to  the 
AS  CVn  binaries  in  Hall's  (1976)  discussion.  It  is  not  a “classical" 

RS  CVn  binary,  so  the  inporunce  of  the  above  paper  Is  clear  only  when 
subsequent  soft  X-ray  data  are  reviewed.  The  soft  X-ray  spectrun  of  the 
“classical"  RS  CVn  binary  UX  Ari  was  fit  with  a plasma  of  10^  “k.  but  it 
could  not  be  fit  with  a plasma  of  solar  abundance  (Walter.  Charles,  and 
aowyer  1979a).  This  gives  additional  evidence  for  an  unOerabundance  of 
heavy  metals  in  RS  CVn  systems.  This  result  was  confirmed  in  two  other 
■classicel"  systems.  RS  CVn  Itself  end  HR  1099,  by  Halter,  Charles,  end 
Sowyer  (19780).  Tne  spectra  of  these  two  stirs  could  not  be  fit  as  nicely 
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as  was  tlie  case  for  UX  An',  but  the  authors  found  that  the  results  were 
consistent  irith  a 10^  thermal  spectra  in  both  cases, 

Walter  et  al.  Cl973a)  felt  that  the  soft  X~ra^  spectra  of  RS  CVn 
binaries  was  consistent  with  an  active  spotted  region.  In  this  model 
they  saw  the  active  region  as  associated  with  flare-like  activity,  which 
1n  turn  would  continuously  supply  suffidient  energy  to  power  a 10^ 
corona.  At  this  temperature  they  computed  the  velocity  Of  the  ions 
to  be  sufficient  to  rapidly  deplete  the  matter.  However,  they  noted 
that  a magnetic  field  of  only  50  Gauss  would  be  enough  to  bottle  up  the 

In  addition  to  the  above,  Walter  et  al.  [1978a]  noted  three  inter- 
esting facts  about  the  new  information  obtained  from  the  X-ray  detectors. 
First,  they  noted  that  1f  UX  Ar1  is  anything  like  a typical  RS  CVn  binary, 
then  from  the  space  density  of  'vlfl*®  pc*^  the  RS  CVn  binaries  contribute 
a full  m of  the  soft  X-ray  background  at  low  galactic  latitudes.  Second, 
they  computed  the  mass  loss  rate  from  their  emission  measure.  The  value 
they  obtained  [10  M_/year)  is  in  reasonable  agreement  with  that  pro- 
posed by  Ulrich  and  Popper  [1974).  Third,  they  reported  that  White, 
Sanford,  and  Weller  [1976)  have  detected  e flare  in  the  X-ray  spectrum 
of  HR  1099  which  was  coincident  with  a radio  flare. 

The  evolutiorarv  status 

The  evolutionary  status  of  the  RS  CVn  binaries  has  been  an 
especially  difficult  problem.  The  mass  of  the  components  is  generally 
cf  the  order  0.5  to  1.5  fi..  This  is  for  both  stars  and  therefore,  with 
few  exceptions,  the  mass  ratio  1s  very  close  to  unity-  The  picture  that 
has  developed  over  Che  lest  10  or  15  years  regarding  binary  star  evolu- 
shown  that  mass  ratios  near  unity  should  be  the 


exception.  RS  CVn  binaries  are  much  too 
case.  Thus  a problem  appears^  either  th 
Hrong  or  the  RS  CVn  systems  do  not  exist 
"bUch-and-nhite"  a statement.  Nore 


nimerous  to  be  an  exceptional 
e binary  evolutionary  scheme  is 
. Obviously,  this  is  nuch  too 
Che  point,  It  would  be  expected 


that  the  theory  needs  seme  adjustment  in  order  for  the  large  nmtber  of 
systems  with  unity  mass  ratio  to  be  explained.  It  Is  important  to  real- 
ize that  the  binary  evolutionary  theory  referred  to  above  explains  the 
evolution  after  the  existence  of  the  binary  pair  on  the  ruin  sequence. 
The  evolution  prior  to  this  point  is  not  well  understood,  and  could  very 
easily  be  the  key  to  the  unity  mass-ratio  puzzle. 

The  similarity  of  the  characteristics  of  the  cooler  component  of 
RS  CVn  to  those  of  T Tauri  stars  led  Hall  {1972)  to  conclude  that  RS  CVn 
was  in  pre-rnaln-sequence  evolution.  This  was  in  contradiction  to  the 
earlier  investigation  by  Field  (19S9),  who  founC  Chat  RS  CVn  as  well  as 
AR  Lac  were  not  In  pre-rnaln-sequence  contraction.  The  age  of  Che  RS  CVn 
binaries  as  Cetermined  by  Hontle  (19731  supported  the  ore-main-sequence 
evolutionary  picture. 


Eventually,  a clearer  case  arose  for  the  post-main-sequence 
evolutionary  status  for  the  RS  CVn  binaries.  This  process  started  with 
□liver's  (1974)  evaluation  of  the  evolutionary  status  and  his  suggestion 
of  a slow  HB55  transfer  by  means  of  a strong  stellar  wind.  Ulrich  and 
Popper  (1974)  proposed  that  by  allowing  for  such  a stellar  wind  it  would 
be  possible  for  the  evolutlonery  status  of  these  stars  to  be  explained  on 
the  basis  of  nonnal  single  star  iDPdels. 

The  pre-main-secuence  status  for  binaries  of  the  RS  CVn  type  was 
almost  conclusively  eliminated  as  a possibility  by  Hall  (1975b),  and 
Slermann  and  Hall  (1976).  Hall  (1975b)  found  that  the  visual  cemipanion  of 


tTie  ecUpsfng  RS  CVn  btnary  WW  Dra  was  a normal  F6  V star.  The  most 
Important  thing  aOoot  1s  was  the  estifliate  of  its  mass.  It  was  less 
massivo  than  either  star  in  the  WW  Dra  system.  This  clearly  would  indi* 
cate  that  the  stars  in  WW  Dra  could  not  Oe  in  ore>main<seouenCB  contrac- 
tion. Hall  also  shoued  that  there  was  sufficient  evidence  of  binary 
(notion  for  the  possibility  of  an  optical  pair  to  be  insignificant- 

Biermanrand  Hall  (1915)  gave  very  convincing  argunents  to  exclude 
the  pre-main-sequence  evolutionary  status.  They  concluded  by  a process 
of  elimination  that  the  BS  CVn  binaries  were  the  result  of  the  fission 
of  a single  star  as  it  leaves  the  main-seouence. 

At  the  time  of  this  writing,  the  consensus  is  that  the  evolu- 
tionary picture  as  presented  by  Popper  and  Ulrich  (1977)  is  the  closest 
to  reality.  They  presented  the  following  circumstantial  evidence  of 
post-maln-seguence  evolution  of  115  CVn  binaries: 

1.  They  are  not  associated  with  regions  of  known  star  formation 

2.  For  the  mass  and  radius  range  the  post-main-seQuence  life-tine  is 
100  times  longer  than  the  pre-main-sequence  life-tio« 

3.  The  larger  radius  is  associated  with  the  more  massive  cooiponent 
A,  WW  Ora  and  HR  1099  have  dwarf  companions  of  lower  luminosity 

In  addition,  they  recalculated  the  age  of  the  RS  CVn  binaries  using 
Nontie's  (1973)  value  for  the  velocity  dispersion  and  Wielen's  (19741 
calibration.  The  results  indicated  that  the  RS  CVn  binaries  are  about 
3x10’  years  old  (i-e.  about  the  same  age  as  similar  mein-seauence  stars). 
In  fact,  they  felt  that  the  age  is  very  slightly  greater  than  a main- 
sequence  star,  and  that  the  RS  CVn  binary  phenomenoidevelops  when  stars 
of  this  general  nature  reach  the  Hertdsprwg  gap.  In  some  cases  they 
proposed  chat  it  would  be  necessary  to  modify  simple  single  star  evplu- 
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The  Intent  of  this  dissertation  can  I 

colors,  separated  in  time  Dy  approximately  o: 
for  each  must  be  obtained  within  a three  month  term  in  order  to  minimUe 
Che  effect  of  a migrating  distortion  wave; 

2.  To  determnne  the  elements  of  the  eclipsing  system  fon  RS  CVn  by  a 
new  technique  which  eliminates  the  problees  introduced  by  tne  oholometric 
distortion  wave. 

[n  meeting  these  two  goals,  several  others  vrill,  by  necessity. 


3.  The  constancy  of  Che  comperlson  star  BIH3S®  2420  will  be  confirmed. 

4.  The  distortion  wave  will  be  isolated  so  that  1t  nay  be  studied  in  a 
future  research  effort. 

5.  The  radial  pulsation  model  for  the  pnotoratrlc  distortion  wave  will 


CrtAPTEH  II 
INSTWWNTATIM 


Tire  Basic  Equipment 

The  Instrunentatlsn  used  1n  collecting  the  data  gresented  in  tMs 
dissertatrlon  may  appear  to  acne  In  the  avant-garde  of  Instrumental 
development  as  antiquated.  This  sentiment  1s  partially  due  to  the 
enonsous  technological  development  during  the  past  decade.  In  addition, 
this  attitude  reflects  a prevailing  view  that  innovation  Is  superior. 
Snobbery  and  high  technological  advance  aside,  the  Important  point  Is 
that  the  data  presented  here  of  US  CVn  are  of  the  highest  quality. 

The  46  cm  telescope  used  In  this  Investigation  1s  located  at 
Rosemary  Kill  Observatory,  6.5  bti  from  Bronson,  Florida.  The  optics, 
telescope  tube,  equatorial  fork  mount,  and  pier  were  built  by  R.  E. 
Brandt.  The  electrical  drive,  slow  motion  controls,  secondary  mirror 


mount,  finder  scope,  counter  weight,  and  setting  circles  were  Installed 
by  the  staff  at  the  University  of  Florida  (primarily,  John  P.  Oliver, 

E.  Whit  Ludington,  and  Eli  SravesI,  The  optical  system  is  an  f/10.5 
Ritchey-Chrdtlen  design.  The  mirrors  of  the  system  are  aluminiaed  and 
Beral  overccated. 

The  phdtcmeter  has  been  described  by  Chen  and  Rehenthaler  (1966), 
Four  dlaphra^  are  avallablai  they  are  small  circular  openings  of 
O.SQO,  Q.2SOO,  0.079  and  0. 039-Inch  diameters  an  a rotatable  disk. 

The  smallest  opening  was  used  in  all  of  the  photometry  presented  here; 
this  corresponds  to  a diameter  of  42. 6 arc  seconds  on  We  sky. 


There  are  three  filters  which  the  user  may  select.  The  light 
passes  through  only  one  filter  at  a tlse,  and  the  observer  loiiet  mamjal- 
1y  move  the  slide  on  which  the  filters  are  mjhted  to  the  desired  filter 
before  a reading  is  taken.  The  ultraviolet  filter  (u)  Is  made  from 
Coming  S6S3,  the  blue  filter  (bj  Is  made  frcm  Corning  5030  and  Schott 
GG13  In  coadilnation,  end  the  yellow  (or  visual)  filter  (v)  Is  made  from 
Coming  3305.  These  filters.  In  union  with  the  same  InstrLiaentatlon  used 
ln  this  photometry  of  RS  CVn,  have  been  found  by  Mirkworth  (1977)  to 
match  closely  the  standard  UBV  systan.  The  photomultiplier  tube  used  in 
this  work  was  an  EHI  9761,  operated  at  900  volts. 

A DC  electrometer  amplifier  designed  by  Oliver  (1976)  was  used 
to  anplliy  the  current  from  the  photomultiplier  tube.  This  amplifier, 
called  PA>10,  has  many  different  gains  which  laay  be  selected  by  rotary 
switches.  This  allows  an  amplification  from  0 to  10^.  The  amplifier 
was  used  with  a time  constant  of  <s  second. 

All  of  the  data  were  recorded  on  a Heath  chart  recorder  with  a 
fiber-tip  pen  opereted  by  e servo-type  motor.  The  observer  may  select 
from  a wide  range  of  speeds  on  this  chart  recorder.  In  all  the  record- 
ings of  BS  CVn,  speeds  of  either  1 1nch/m1n  or  0.5  1nch/m1n  were  used. 

from  a fligltal  integrator,  'nils  instrument,  referred  to  as  DA-10,  was 
designed  and  built  by  John  P,  Oliver.  It  Integrates  the  signal  by 
counting  the  pulses  from  a vol  taga-to-frequency  converter  1n  the  PA-10 
amplifier.  The  observer  has  the  option  of  irany  integration  times  from 
1 second  to  100  seconds.  The  Integrations  are  started  manually, 
but  are  stopped  and  displayed  by  the  Internal  electronics. 


Fig.  2.'Sc1ifijiiat1c  diagram  of  the  instrocnentatlori  used  to  obtain 
the  ohotonetrlc  data  of  RS  CVn  at  Rosemary  Hill  Observatory  during 


In  Figures  3 end  4.  the  relative  gains  for  settings  of  the 
fine-gain  and  coarse-gain  switches  are  presented,  respectively.  Fmm 
inspection  of  these  figures,  it  Is  dearly  evident  that  the  short-tenn 
stability  and  the  long-ten  stability  of  the  amplifier  are  very  good. 
In  Fig,  4 a difference  in  the  level  of  the  coarse-gain  is  evident  from 
1975  to  1976.  This  is  of  no  significance  in  differential  photometry 
such  as  that  presented  here  on  US  Canum  Venacicarom, 
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. can  he  achieved  with  inadern  photoelecCric 
• the  earth's  atmosphere  more  than  by  any 
s atmosphere  is  not  completely  transparent; 


it  does  not  even  have  the  same  transparency  at  a11  wavelengths  of 
light-  The  problem  of  measuring  the  intensity  of  stars  Is  further 
complicated  by  the  random  and  unpredictable  fluctuations  in  the  trans- 
parency of  the  atmosphere  as  a function  of  time,  and  the  Sdmewhat  pre- 
dictable change  as  a function  of  aimass. 

There  are  means  of  reducing  the  effects  of  a real  atmosphere 

star's  light.  The  most  important  is  to  measure 


the  light  differentially.  By  this  it  is  meant  that  a star  of  constant 
brightness,  called  the  eonpariaon  star,  is  chosen  as  a reference-  The 
brightness  of  the  program  star,  called  Che  variable  star,  is  measured 
under  the  same  atmospheric  conditions  as  the  brightness  of  the  compari- 
son star-  If  the  atmospheric  conditions  are  exactly  the  sans  when  the 


brightnesses  of  the  two  stars  are  measured,  and  the  comparison  star  is 
Indeed  of  exactly  constant  brightness,  then  any  variation  in  the  ratio 
of  the  varlaOle  star  brightness  to  the  comparison  star  brightness  Is 
due  solely  to  the  variable  star. 

Of  course,  it  Is  not  possible  in  practice  to  guarantee  that  the 
conditions  are  exactly  the  same  during  the  measurement  of 
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conpar'lson  star  brightness.  In  mklng  these  calculations  it  will  be 
possible  to  correct,  in  part,  for  the  atmospheric  extinction  which  has 
been  ntnimizeb,  but  not  eliminated  by  the  proper  choice  of  a comparison 

It  was  pointeO  out  in  Cnapter  II  tnat  the  majority  of  the  data 
exists  in  two  forms:  chart  records  and  tabulations  manually  recorded 
from  the  DA- 10  digital  integrator.  The  chart  records  maoe  prior  to  S, 

6 March  1976  were  read  and  recorded  in  a tabulation  form  exactly  like 
that  used  on  and  after  S,  6 March  1976.  The  keypunch  shop  of  the 
northeast  Regional  Data  Center  was  engaged  to  punch  these  data  on  stan- 
dard computer  cards,  in  a format  called  the  "Ludlngton-Oliver'  format 
{see  Table  4).  A computer  program  called  RHRED  (PHotoelectric 
REtkjctiohs)  was  written  to  accept  this  format  for  the  raw  data  and  to 
do  the  reduction  to  the  instrumental  system  (except  for  the  inclusion 
of  the  effects  of  atmospheric  extinction).  The  application  of  the  ef- 
fects of  a real  atmsphere  to  the  partially  reduced  output  of  PHRED  is 
handled  by  a computer  program  called  DEXTOR. 

Introduction 

PKRCD  is  a computer  program  that  is  written  in  the  PLl  program- 
ming language.  It  is  a very  versatile  program  for  reducing  photoelec- 
tric observations  to  the  inetrumental  system.  It  has  far  vaster  capa- 
bilities and  fewer  limitations  tnan  any  other  available  reduction  pro- 
gram. It  is  limited  by  the  size  of  the  computer  used  rather  tnan  by 
the  observations.  It  oan,  for  example,  handle  the  observations  of  over 
600  different  objects  in  JOO  different  filters  taken  during  one  night 
of  observing.  PHRED  can  be  instructed  tc  produce  a plot  for  any  object- 


evaluating  tne  quality  of  tne  observatians.  In  addition,  tiie  uaer  can 
instruct  l>HfE0  to  calculate  the  principal  entinction  coefficients,  the 
auts1de>the>aCnosphere  inagnicude  and  the  standard  errors  of  each. 
Furthermore,  the  program  is  not  restricted  to  variable  star  data.  It 
can  also  do  the  reduction  of  standard  stars  to  the  instrumental  system. 
PHREO  does  a11  rtortnal  reduction,  except  atmospneric  extinction  and  the 
reduction  to  the  standard  system.  For  example,  it  corrects  for  a zero 
offset  in  the  data>gathering  system.  It  converts  to  heliocentric  Modi- 
fied Julian  Date.  It  applies  gains,  unless  the  user  is  using  a pulse 
counting  system,  in  which  case  the  gains  will  be  ignored.  All  in  all, 
the  user  has  a great  deal  of  control  over  the  operation  of  the  program, 
and  very  few  assumptions  have  been  made  about  the  way  the  data  must  be 
taken,  or  the  organization  they  must  be  in  for  reduction  by  this  pro- 
gram. Of  equal  importance  to  the  user  of  PHRED  is  the  detailed  log  of 
the  reduction  and  Che  warning  messages  to  help  the  user  gain  the  de- 
sired results. 

Prodram  Control 


The  user  instructs  the  program  to  do  certain  things  by  use  of 
Control  Cards.  These  are  records  that  start  with  an  asterisk  (*),  The 
asterisk  is  followed  by  one  letter  of  the  alphabet.  This  letter  deter- 
mines the  option  that  is  to  be  used  or  the  type  of  parameter  cards  Chat 
will  follow  it-  Any  characters  after  the  first  letter  will  not  be  re- 
cognized by  the  program. 

The  first  step  in  using  PHREO  1s  to  have  the  data  prepared  on 
cards  1n  the  Ludington-OUver  format.  The  outline  of  this  format  is 
given  in  Table  4,  Any  item  on  which  a deflection  can  be  taken  is  de- 
fined by  the  user  with  an  object  code.  This  includes  sky  deflections. 
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d&rK  current,  zero  offset,  var'fable  star,  standard  star,  etc.  Hiere 
strictions  on  mMcN  of  the  ZSS  characters  of  the  EBCDIC  character  set 


on  a card  for  three  filters  to  be  coded.  If  more  are  required,  the 


data  are  continued  on  the  next  card  by  copying  the  object  code.  The 
gain  Is  an  instrumental  setting  Chat  causes  the  Input  from  the  photo- 
muUlpler  tube  to  be  amplified  by  a specific  amount.  This  must  be  ac- 
counted for  1n  order  to  arrive  at  the  true  value  before  amplification. 
The  code  that  Is  used  on  the  Lud1nglon-01 Iver  format  is  comprised  of 


while  the  second  is  for  a setting  of  a f1ne-ga1n  stvitc 
not  make  any  distinction  between  the  two.  It  reads  th 
as  a single  code;  therefore,  the  user  must  supply  (In 
section}  the  value  for  this  combination  of  coarse  and 


settings.  The  deflection  thi 
amplification  indicated  by  tl 
to  five  digits.  This  should 
photoelectric  photometry.  Tf 
bedded  blanks.  A deflection 
from  the  photoelectric  devici 
offset  and  must  be 
are  to  be  achieved 


iree  characters 
control  card 
i-ga1n  control 

Is  recorded  will  of  course  Include  Che 
gain  code.  The  deflection  Is  restricted 
1 sufficient  for  the  most  majority  of  the 
deflection  fields  must  not  Include  em- 
i possible  even  when  the  signal  current 

amoved  from  other  deflections  1f  the  proper  results 
Therefore,  PHREO  will  recognize  an  object  code  'Zb' 


as  the  zero  reading,  It  will  ignore  Ih 
average  of  all  such  readings  for  that  n 
from  all  deflections  before  the  gains  a 


gain  satting  used  , 
gnt  and  substract  tf 
e applied.  The  usei 
ir  each  obsarvatlon.  1 


Julian  Date  (irtilcti  was  definea  by  tba  lAb  in  1973)  is  recorded  on  each 
card;  this  is  actually  the  interjer  oart  of  the  iloOifled  Julian  Date. 

The  Universal  Tima  (UT)  is  also  recorded  for  each  deflection.  Sij 
colunirs  are  used  for  each  UT.  The  first  two  digits  are  for  the  hours, 
the  second  two  digits  are  for  the  minutes,  and  the  last  two  digits  are 

PHffiD  has  been  designed  to  perform  several  tasks  for  the  user. 
PHRED  can  process  one  night  or  mitiple  nights  of  observations  in  a 
single  submission.  The  first  night  of  observations  begins  with  a set 
of  control  cards;  these  might  Include  the  'LiBRART,  *FILTERS,  *CODES, 
»S0RT,  'VARIABLE,  *i>L(JTS,  'BEGIN,  and  'TITLE  control  cards.  The  last 
two.  however,  have  a special  importance.  The  ”9ESIN  must  be  the  first 
card  in  the  control  section  and  the  *nTL£  rajst  be  the  last  control 
card  before  the  data.  The  other  control  cards  may  be  in  any  order. 

There  are  three  additional  control  cards  which  are  of  special  importance. 
These  are  the  'DATA,  'ADD,  and  'END  control  cards.  'DATA  must  follow 
the  title  cards  and  proceed  the  data.  wENO  must  follow  the  data.  'ADD 

than  to  be  considered  as  the  original  data.  For  ejample,  the  observer 
may  have  neglected  to  take  a comparison  star  deflection  at  the  very  be- 
ginning of  the  night.  While  this  will  not  prevent  the  normal  operation 
of  PHRED,  it  will  require  extrapolation.  To  avoid  this,  the  user  can 
add  an  artificial  deflection  with  the  'ADD  control  cord.  After  the 
original  set  of  data  has  been  ended  with  an  *END  card,  the  user  inserts 
the  'ADO  card  followed  by  as  many  cards  of  additional  data  (in  the 
Ludington-Oliver  format)  as  he  desires.  The  additional  data  are  also 
ended  with  an  wENO  card.  The  next  night  of  observations  will  then 


follow,  beginning  with  the  •BESIN  CARO,  It  Is  not  necessarji  to  re- 
peat all  control  cards  for  the  second  night  it  they  were  already  used 
in  the  first  night's  control  card  section.  That  is,  the  second  night 
will  be  processed  with  the  sane  instructions  as  for  the  first  night, 
unless  the  user  changes  then.  For  example,  it  Is  unlikely  that  the 
gains  will  have  exactly  the  same  xalues  from  night  to  night.  Therefore, 
the  'GAIN  control  card  will  probably  be  necessary  in  the  second  night's 
control  section.  On  the  other  hand,  the  stars  that  were  observed  and 
the  filters  that  were  used  will  probably  be  the  sane;  therefore,  the 
user  does  not  need  the  •LI9RART,  'COOES  and  'FILTER  control  cards  in 
the  control  card  section  for  the  second  night.  Four  of  the  control 
cards  work  in  a "flip/flop"  fashion,  i.e,,  they  are  "off"  until  they 

and  then  they  are  "off,  etc.  These  are  'ItSGNITUOES,  'INTERPOLATE  ALL 
SO,  *S0RT,  and  'VARIABLE. 

In  Figure  5 is  a listing  of  the  JCL  cards  for  using  PNRED. 
Notice  that  this  assumes  the  existence  of  a load  nodule  for  FHRED. 


Control  Cards 

Listed  hare  are  the  control  cards  and  their  meanings: 

'ADO  --  Adds  additional  Data  Cards  to  those  already  read  and  writes 
then  into  the  Raw  Data  file  (RAVJD}.  The  additional  data  shoulO  be  ended 


witli  *END.  The  Lutflngton-Ollver  format  ia  used  far  additianal  data. 

night  of  data.  Each  time  this  card  Is  encountered,  the  oreyious  night 
is  known  by  PHRED  to  be  ccsnplete  and  is  then  processed. 

•BOTTOM  - This  card  should  be  the  last  card  in  the  input  deck. 

It  will  insure  the  processing  of  the  last  night  of  data. 

•COOES  FOR  OBJECTS  — The  following  parameter  cards  define  the 
object  codes.  A blank  code  terminates  object  code  input.  The  codes  are 
two  characters  long,  left  justified  in  a fcur>co1miin  field.  For  eiample, 
V SV  C SC  K SK 

is  the  default  if  no  codes  are  defined.  Any  code  used  in  the  data  but 
not  defined  will  cause  a warning  message  to  be  printed  in  the  LOG  file, 
and  the  data  will  be  ignored.  Any  code  that  is  defined,  but  for  which 

stellar  cades  should  be  odd  (1st,  3rd,  5th,  etc.),  sky  codes  should  be 

•VARIABLE  option  is  used. 

•DATA  --  This  card  signals  the  beginning  of  the  data  card  ob- 


ten  into  the  Raw  Data  file  (KAWO).  The  •DATA  card  must  follow  the  last 
title  card  (see  ‘TITLE). 

•END  “ This  card  is  used  to  end  a group  of  data  cards.  Either 
the  original  data  or  any  and  all  additions  to  it  are  ended  with  ‘END. 

•FILTER  CODES  - The  following  parameter  cards  define  the  filUr 
codes.  A blank  code  terminates  the  set.  The  codes  are  one  character  In 
length,  left  justified  in  two-column  fields.  For  example, 

is  Che  default  if  no  codes  are  defined. 


*GAIKS  The  fo11owin9  paraineter  cerds  define  the  gain  codes 
and  their  values.  A blank  code  and  value  terminates  the  set.  The 
codes  are  three  characters  In  length,  uh11e  the  value  can  be  as  much  as 
12  columns  long.  The  code  and  its  value  are  separated  by  one  space  (an 
epual  sign  or  any  other  single  character  may  be  used).  Therefore,  to- 
gether they  comprise  16  colums-  For  ekample. 

D04vO.SO«9  E02*2. 501234567  A06=  -0.49S9 

The  default  for  this  control  card  is  to  apply  a gain  of  zero  magnitudes 
to  all  the  data  points.  A message  will  be  written  into  the  LOG  to  indi- 

rwARDCOFY  — This  card  is  included  if  a punched  deck  (or  other 
machine  readable  copy)  is  desired  of  the  reduced  data.  The  parameter 
card  can  be  used  to  control  the  contents  of  this  output  file.  The 
parameter  card  format  is  a six-character  field  in  columns  one  through 
six.  There  are,  at  the  present,  only  three  possibilities  for  the  user 


'ALL"  will  cause  the  total  of  the  iiSXY  printed  file  to  be  put  into  the 
MCK  output  file.  "HOME"  will  cause  nothing  to  be  put  into  the  DECK 
output  file,  while  "V/C"  will  cause  only  the  variable  minus  conrarison 
data  to  be  put  into  the  DECK  output  file, 

'IKTERPQUTE  ALL  SKY  — IKually  an  observer  will  take  deflections 
on  a star  and  on  the  nearby  sky.  This  will  allow  him  to  subtract  the 
sky  from  the  star  deflection  and  be  left  with  Just  the  deflection  due 


code  Interpolation  of  sky  deflections.  TMs  means  that  If  there  are 
sky  observations  only  with  the  companison  and  none  for  the  variable, 
the  user  must  Include  this  card.  Without  this  card,  the  projram  sub- 
tracts the  interpolated  value  of  an  even  code  from  the  value  from  the 

b1ned  and  copied  in  place  of  the  original  even  code  values. 

*JD  OVERRIDE  — The  following  parameter  card  defines  a five- 
digit  number  that  will  override  the  Integer  Wodifled  Oulian  Date  on 
the  data  cards.  It  must  not  have  a decimal  point  or  sign.  For  example. 


42451 

Is  a valid  parameter  card.  The  default  Is  to  use  the  integer  Modified 
Oulian  Cate  on  the  data  cards.  Note  that  the  Integer  Modified  Julian 
Date  coded  on  each  card  Is  used  to  conpute  the  time  for  each  of  the 
deflections  on  the  card.  Therefore,  the  use  of  this  control  card  op- 
tion will  cause  all  the  cards  to  have  the  same  Julian  Date,  even  if  the 
originel  data  cards  for  a single  night  contained  more  than  one  Julian 
Date.  The  default  is  reestablished  at  the  beginning  of  a new  night  of  data. 

n-EniNa:0H  COEFF  - if  the  user  wishes  to  have  extinction 
coefficients  calculated,  he  may  specify  on  the  following  parameter 
card(s)  which  star  (or  sky)  and  f11ter(s)  he  wishes  used.  There  are 
no  limits  to  this  oomOinatlon  of  object  code  and  filter  code.  The  for- 
mat for  the  parameter  card  Is  a four-column  field  with  the  object  code 
in  the  first  two  columns  and  the  filter  code  in  the  third  column  {the 
last  column  Is  blank).  For  example. 


vjsv  5VV  cev  SCV 

will  produce  extinction  coefficients  for  the  variable  star,  for  the 
sky  readings  next  to  the  variable,  for  the  ccmoarisan,  and  for  the 


coniparlson  sky,  all  for 


tke  visual  filter  (assuming  tlie  object  and 
filter  codes  defined  In  ttie  examples  for  rcODES  FOR  OBJECTS  and 
FILTIRS}.  The  default  is  to  not  calculate  any  coefficients. 

•LIBRARY  STARS  --  The  following  parameter  cards  contain  the 
"vital  statistics"  on  the  stars  being  reduced.  The  format  for  these 


cards  appears  below: 


B0-15“l734 


RA 


EXAMPLE  eommis 

C»  01-02 

BD-15  1734  03-12 

7.2S33333333  14-2B 
-15.S41imn  29-43 


FOWAT 

A{2) 

A(10) 

F{15,10) 

F(I5,I0) 


Rote  that  when  punched  on  the  cards,  the  right  ascensions  and  declina- 
tions must  be  converted  to  hours  and  decimal  hours  and  degrees  and  deoi- 


mal  degrees.  The  last  card  must  be  a blank  card  to  terminate  the  set. 
If  no  library  parameter  cards  are  used,  then  the  •X-EniNCTION  COEFF 
control  card  option  should  not  be  used. 


•HASHITUDE  - If  this  card  is  used,  the  output  will  be  in 
magnitudes  instead  of  intensity  units.  The  default  is  output  in  in- 
tensity units. 

*N  FILL  SCALE  — The  following  parameter  card  will  give  the 
fuHscele  value  expected  for  the  data.  This  value  will  then  be  used  to 
normalize  all  the  deflections.  This  means  that  before  anything  else, 
the  values  punched  on  the  data  cards  will  be  divided  by  this  number. 

The  number  is  read  in  from  the  first  six  columns  of  the  parameter  card 


using  an  F(6)  formiat.  For  example. 


are  a11  va1td  fu11~sca1e  values.  The  default  va'iue  fs  IDOOQ. 

voeSERVATORV  — The  foUowing  parameter  card  is  used  to  read 
in  the  longitudei  latitude,  and  (optionally)  the  station  code  for  the 
oOaervatory  at  which  the  data  were  taken-  For  example, 

♦H. 5866666666  f29. 377777777  WO 

is  the  default  and/or  the  parameter  card  for  Rosemary  N1I1  Observatory, 
It  Is  read  by  (F(15,10).F(15,10),X(5),A{3))  PLl  foniiat. 

Longitude  is  a number  between  >180,0  and  fl60.0, 

*PLOTS  “ The  following  parameter  oards  are  used  exaotly  like 
the  one  for  wK-EXTinaiOh  COEFF.  They  specify  which  object-filter 
combinations  the  user  wants  plotted.  The  codes  and  formats  are  the 
same  as  for  the  *)!-EXTINCnON  COEFF  parameter  cards.  The  plots  are  in 
magnitudes  as  a function  of  local  time.  No  plot  is  the  default. 

•SORT  ” This  card,  when  used,  will  cause  the  sorting  of  every 
object-filter  combinetion  into  time  order.  Sorted  data  are  necessary 
for  proper  functioning  of  PHRED. 

•TITIE  — This  card  signals  the  beginning  of  the  title  cards. 
The  first  of  these  cards  will  be  used  as  a heading  on  all  pages  per- 
taining to  this  night  of  data.  The  rest  (and  there  is  no  limit}  will 
be  printed  an  the  heading  of  the  RAWD  print  file.  There  must  be  a 
•TITLE  card  and  at  least  one  card  following  it.  the  •DATA  card  must 
then  follow  the  last  title  card. 

•VARIABLE  — If  the  data  are  arranged  for  it,  the  user  can 
have  the  data  treated  as  "variable  star  data."  In  this  case  the  com- 
parison star  will  be  "removed"  from  the  variable  (to  get  relative 
brightness).  This  is  done  by  dividing  the  variable  star's  reduced  in- 
tensity by  the  Interpolated  conparison  star's  reduced  intensity.  If 
•MAGNITUDE  has  been  used,  the  process  will  be  the  difference  of  the 


reduced  Tagnltudss.  Th«  object  codes  loust  be  listed  after  the  'CODES 
FOR  OBJECTS  control  card  In  this  order;  variable  star,  sky  variable, 
comparison  star,  sky  con^arlSDn.  another  star,  another  sky,  etc. 

■*x  TIME  OFFSET  — The  following  parameter  card  can  have  op  to 
a lS-d1g1t  number  {read  by  F(15,14}).  This  number  Is  added  to  the 
Julian  data  computed  for  each  deflection.  Its  units  must  therefore  be 

DEXTOB 

Prooram  Outline 

DEXTDR  Is  a program  written  In  PLl  which  applies  differential 
extinction  coefficients  to  a dataset  like  that  produced  by  PNRED. 

If  the  Instrumental  system  is  nearly  the  same  as  Che  standard 
UBV  system,  a simple  linear  set  of  equations  can  accurately  describe 
the  effects  of  a real  atmosphere.  The  equations  normally  used  are 

tla)  v=  . v-k-X, 

(lb)  (b-v)®  . ((b-v)-k,-X)/Cl*k|-X),  and 

(lc)  (u-b)°  = ((u-bj-kj.Xj/dtk^.X), 

where  the  superscript  "o"  Indicates  the  outsice-the-abnosphere  color 
or  magnitude  of  the  star,  X is  the  airmass,  v.b,  and  u are  the  nagn1< 
tudes  of  the  star  as  viewed  through  the  visual,  blue,  and  ultraviolet 
filters,  respectlyely,  and  k,  k^,  kj.  kj,  and  kj  are  the  extinction 
coefficients.  For  a Cetailed  development  of  the  above  equations  see 
Kardie  (1962). 

The  above  equations  apply  to  the  measurement  of  one  star,  but 
the  data  can  better  be  reduced  differentially,  rather  than  as  in- 
dividual magnitudes  of  the  variable  and  comparison  stirs.  In  order  to 


differential  oktlhction  computations.  Start  by  defining  a difference 


operator  D(..0  which  is  equal  to  the  variable  minus  the  comparison 
in  magnitudes.  The  above  equations  can  Chen  be  defined  for  the  vari- 
able star  and  for  the  comparison  star,  it  can  be  shown  that  equations 


(a)  0(v)“  » Il(Y)-k-D(!(), 

(2B)  D(b.v)»  . °(b-T)-fci-gU)-4-n(f(b-v))_ 

(20)  D(u-bl°  = 0(»b)-k;-l)(X)-4-D(X-(u.b))  ^ 

(l*k^-X,)  (lrkj-*|.l 

where  Che  symbols  have  the  same  meanings  as  for  equations  tl),  with 
the  adoitiohs  of  X equals  the  airmass  of  the  variable  star,  and  X 
equals  the  airmass  of  the  conparison  star. 

The  equations  (2)  above  are  the  ones  which  are  used  by  D2XT0R 
in  computing  Che  instrumental  differential  magnitudes  of  the  variable 
star.  The  eatinccion  coeffioients  must  be  ccmpuced  separately  and 
supplied  to  the  DEXTOR  program  via  control  cards. 

In  the  domputations  perfprmed  Py  DEXTQR,  a simplifying  assump- 
tipn  has  been  made  about  the  values  of  0(X-(b-v))  and  0(X-(u-b)].  It  is 
assumed  that  the  differenoes  can  be  replaced  by  the  appraximations, 

13a)  D(X-(b-,))v\^»d,p(t,.,)_ 


The  error  introduced  by  this  approximation  is  small  if  D{X)  Is 
small.  For  D(X)  ■ 0.01  the  error  ih  D(B-v)®  is  only  0^04  for  typical 
extinction  coefficients,  ho  error  was  introduced  in  the  reduction  of  the 


adopted  for  k^- 


Extinction  Coefficients 


An  aversge  value  for  each  of  the  extinction  coefficients  was 
determined  Oy  the  method  of  weighted  least  squares.  The  data  used 
were  those  from  the  comparison  and  check  stars.  The  evaluation  of  the 
coefficients  was  based  on  equations  similar  to  equation  (la).  The 
evaluation  of  these  primary  ertinction  eoefficlehts  was  perfoniied  in 
the  PHRED  computer  program  for  each  night  of  data.  The  results  of 
PKREO's  computations  were  then  weighted  according  to  the  precision  of 
PHREO's  least  squares  determination,  and  the  final  extinction  coef- 
ficients were  obtained  by  an  additional  application  of  the  method  of 
least  squares.  The  values  adopted  for  application  by  DEXTdR  to  all  the 
data  (I97a  and  1976)  were: 


k,  • 0.667,  and 


The  values  for  kj  and  k^  were  adopted  based  on  the  closeness  of  the 
computed  values  to  the  expected  values  given  by  Hardie  (196E). 

The  details  of  the  computations  are  not  critical  to  the  quality 
of  the  results.  The  effects  of  differential  extinction  are  small. 

In  feet,  the  average  value  of  the  correction  on  a typical  night  was 
0.001  magnitudes  In  the  D(v)^  value. 


After  the  data  were  punched  onto  computer  cards,  the  cards  were 
PHRED,  It  was  necessary  to  first  punch  all 


collected  for  processing  by 


of  the  control  cartfs  for  PHfiED  and  yeHfj  that  all  the  supplemental 
data  transferred  to  the  program  via  the  control  section  were  correct 


Mhen  the  above  conditions  were  satisfied,  the  PHUED  portion 
of  the  reduction  was  complete.  The  computer  listings  are  not  included 
here  because  they  are  bulky  and  do  not  easily  conform  to  the  reguired 
format.  Copies  are.  however,  available  to  interested  investigators. 
Correction  of  Errors 

Since  there  were  a large  number  of  data,  it  was  impossible  to 
verify  all  nioidrers  in  the  tabulations,  reduction  and  pnintout.  However, 
It  wes  important  to  investigate  all  significant  human  errors  so  that 
they  would  not  propagate  through  the  remelhing  reduction. 

The  first  step  to  reduce  human  errors  was  taken  during  the 
keypunching.  The  tabulated  data  were  punched  and  then  verified  by 
the  keypunch  staff.  Second,  a progrmn  was  written  to  scan  Che  punched 
cards  and  Co  flag  changes  in  the  deflection  value  Chat  were  obviously 
too  large  or  were  changes  in  time  sequence.  This  would  locate  many 
transpositions  of  digits  in  the  tabulation  process,  or  keypunch  errors. 
The  program  was  only  partially  successful.  Third,  the  reduced  data 
were  converted  by  use  of  the  equation 
(4)  HJO  = «:825.042  * 4.797BS5E, 
to  a common  phase  cycle  and  the  light  curve  plotted.  The  linear 
ephemeris  above  is  due  to  Catelano  and  Rodono  (1874).  Any  datum 
which  was  obviously  discrepant  was  carefully  checked  for  errors  in 
tabulation,  kaypunching,  or  data  reduction.  Fourth,  to  check  the  re- 
duction process  performed  by  PKRED,  a ssn^le  set  of  calculations  was 
perfonned  manually  with  an  HP-45  nand-held  calculator.  The 


sample 


Cilculatlons,  In  every  respect,  sgreeO  nitli  ttie  vslees  fnw  the  PHBED 
program  to  wtthin  Che  expected  accuracy. 

The  majority  of  effort  mss  devoted  to  the  third  method, 
descrihed  above,  for  detecting  errors.  More  than  60  individual 
data  were  carefully  checked  in  detail  [requiring  10  full  days  of  the 
writer's  undivided  attention).  The  majority  of  the  errors  detected 
occurred  at  the  time  the  data  were  tabulated  for  keypunching.  The  cor- 
rections could  easily  be  made  and  the  data  adjusted,  in  a few 
cases  it  was  possible  to  identify  the  source  of  the  error  (miss-identi- 
fication of  star,  clouds,  telescope  drift,  etc.),  but  the  information 
necessary  to  correct  these  errors  was  not  available.  In  these  few 
cases  it  was  necessary  to  remove  the  offending  data  frcm  the  collection 
in  order  to  avoid  contamination  of  the  results  by  these  known  errors. 

In  only  one  case  was  a datum  renmved  without  an  explanation  that  was 
con*1etely  sitisfactory  to  this  writer.  The  first  star  deflections 
(v.b.and  u)  were  all  significantly  too  high  and  the  first  variable  star 
deflections  (v.b  and  u)  were  significantly  low  on  the  night  of  1.  2 
April  1976.  for  the  remainder  of  the  night  the  data  were  very  good, 
therefore  these  first  measurements  were  roBved  from  consideration. 
Check  Star  Data 

The  ratio  of  the  check  star  (B[H35'’24W)  intensity  to  the 
comparison  star  intensity  should  not  be  time  dependent  if  the  stars 
are  both  of  constant  brightness.  This  allows  the  observer  to  have  a 
'check'  on  the  assumption  that  the  brightness  of  the  coinparison  star  is 
constant,  and  therefore  that  the  variable-to-comparison  ratio  is  a true 
measure  of  the  variation  in  the  variable  star's  brightness  only. 

Occasionally  during  the  course  of  the  observing  program. 


TABLE  5 
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plot  of  ctie  Instruirerital  differential  magnitudes  is  presented  as  a 
function  of  Modified  Julian  Date,  Notice  that  the  plot  Is  not  con* 
tinuous  in  time;  this  is  because  the  observing  was  carried  out  only 
during  the  spring  months  of  1975  and  1976,  The  data  plotted  in  Fig.  6 
appear  in  tabular  form  in  Table  S. 

The  data  for  RS  CVn  obtained  with  the  Instrumentation  des- 
cribed in  Chapter  II,  during  the  1975  and  1976  observing  seasons 
from  March  to  Hay.  and  reduced  by  PHRED  and  DEXTOR,  appear  in  tabular 
form  as  Instrumental  differential  magnitudes  in  Che  appendla. 
Transformation  to  Standard  System 

The  goals  of  this  research  do  not  reguire  transformation  to  tne 
standard  system.  Differential  measurements  with  the  same  ccngiarison 
Stan  and  subsequent  use  of  HIHK  in  the  analysis  (see  Chapters  IV  and  V) 
relaies  the  requirements  for  standardisation  1n  the  data.  For  these 
reasons,  only  a small  effort  was  made  to  obtain  transformati on  coeffi- 

The  standard  stars  used  were  observed  on  two  nights  during  the 
caurse  of  the  observing  program.  The  transforaetion  eouatlons  that 
apply  are  given  as  follows: 

(Sa)  6-V  ■ Aj  ♦ A2<b-v)® 

(5b)  U-B  = A3  f Aj(u-b)“ 

(5c)  V-v“  = Aj  * Ag(b-v)® 

The  values  for  the  A's  in  equations  (5)  were  computed  by  the 
method  of  least  squares  using  the  data  obtained  from  the  photometry  of 

not  been  transformed  with  these  coefficients  because  the  accuracy  is  so 


trafisftsmatJon  coefficienu 


> 1.06  £ Q.15  (tlie  estimated  errors  are 


Hare  accurate  transfomatlons  can  be  found  In  Hartwortn’s  {1977} 
dissertation.  His  ohottanetry  used  virtually  Identical  instrumentation. 


CHAPTER  Itf 

DATA  ANALYSIS  PROCEDURE 
Introduction 

Die  analysis  of  Che  data  was  carried  out  on  an  ANCUSHL  470  cdm- 
outer  at  the  Northeast  Regional  Data  Center  located  on  the  Gainesville 
campus  of  the  University  of  Florida.  The  processing  was  accomplished 
ivitn  several  programs  written  in  FORTRAN  IV.  One  program  IHINK).  iHilch 
was  written  by  0.  6.  Hood  (1971,  1972),  oodels  an  eclipsing  binary  system. 
The  remaining  programs  were  written  by  the  author  except  for  the  least 
square  subroutine  (used  1n  the  Fourier  program)  which  was  written  by 
N.  L.  Cohen.  The  details  of  the  operation  will  be  discussed  after  an 
butHne  bf  the  HAVE  procedure  1s  discussed. 

Outline  of  HAVE  Procedure 

A simplified  outline  of  the  procedure  appears  In  F1g.  7.  The 
details  of  these  comoutations  are  elltrlnated  so  tnat  the  process  is 
clearer.  The  procedure  that  this  figure  represents  was  used  as  an 
in>stream  procedure  and  was  called  HAVE.  This  term  will  be  used  In  the 
remainder  of  this  dissertation  to  refer  to  this  seguence  of  computer 
programs  which  was  used  to  snalyae  the  data.  Die  HAVE  procedure  starts 
basically  with  tbe  HINR  program  which  both  Improves  the  eclipsing  binary 
elements  and  produces  e theoretical  light  curve.  The  input  data  set  to 
HINK  (rdiidr  it  attempts  to  solve)  Is  celled  OROQS.  Dris  term  refers  to 
the  Old  Revised  WServatlons  from  a previous  HAVE  run  or  to  observations 


Fig.  7..Qutlir 
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thwiselves  wher  ttie  procedure  Is  Initiated.  In  addition  to  the  OllDBS 
the  user  also  supplies  a set  of  control  parameters  to  UWK.  After  five 
HINI!  iterations  or  after  convergence  (whichever  cones  first)  the  WINK 
program  produces  a dataset  called  CALCM.  This  is  the  theoretical  light 
curve  for  the  eleeients  at  the  conclusion  of  the  solution  portion  of 
S/IKK.  The  term  (CALOi)  refers  to  the  fact  that  it  is  necessary  to  pro- 
duce Che  calculated  light  curve  in  magnitudes  in  order  to  include  Che 
Quadriture  Magnitude  (QH)  parameter  in  the  values  produced.  This  data- 
set and  the  observations  (OBS)  are  used  in  the  ROMC  program  which  is 
next  in  the  seguence.  fiOMC  converts  the  magnitudes  of  the  CALCM  dataset 
to  intensities  and  subtracts  then  from  the  CBS  dataset.  The  result  is 
the  CMC  dataset,  which  1s  a representation  of  the  distortion  wave.  In 
addition,  iiCHC  produces  a CALC  dataset  which  is  the  converted  CALCM 
dataset.  CALC  is  used  later  for  plotting  purposes.  Next  in  the  sequence 
Is  the  program  FDLIRriT.  It  uses  the  CMC  dataset  to  deUraine  the  coef- 
ficients of  a truncated  Fourier  series,  by  tlie  method  of  least  squares, 
which  represents  the  distortion  wove.  The  Fourier  representation  is 
then  subtracted  fnwi  the  DBS  dataset.  The  result  of  the  subtraction  Is 

metric  complication.  This  dataset  is  called  ROBS  (tevised  OBS).  FOUR- 
FIT  does  additional  computations  whicn  are  used  for  plotting  purposes. 

It  computes  a dataset  from  the  truncated  Fourier  series  at  closely 
spaced  phases.  This  dataset  is  called  FOHC.  Lastly,  FOUfiFIT  calculates 
the  difference  between  the  OMC  dataset  and  the  Fourier  representation  of 
it.  This  is  called  TOMC. 

The  last  significant  program  in  the  MAVE  procedure  is  called 
LCFLOT,  This  program  plots  the  datasets  produced  in  the  preylous  steps 


78 

of  the  procedure  for  evaluation  by  the  user.  A sample  of  the  plots 
produced  appears  In  Fig.  8. 

The  VJINA  Program 

This  program  has  beeh  widely  used  to  determine  the  elements  of 
eclipsing  binary  systems.  The  progrmn  was  written  by  0.  B.  Wood  (1971, 

1977}  and  has  been  updated  by  eight  status  reports  (private  comnuhlcations} . 
The  model  Is  that  of  tnaxlal  ellipsoids.  The  parameters  of  the  model 
which  are  pertinent  to  this  Investigation,  their  WINK  codes  and  meanings 


WINK  PARAMBTIItS 


Symbol  rieaning 


orbital  inclination 

e is  the  orbital  eccentricity  and 

u is  Che  longitude  of  perlastron 

time  (or  phase]  of  conjunotion 

linear  limb  darkening  coefficient  for  star  A 

linear  limb  darkening  coefficient  for  star  B 

reflection  albedo  of  star  A 

reflection  albedo  of  star  B 

guadrature  magnitude 

effective  temperature  of  star  A 

effective  temperature  of  star  8 

radius  of  a sphere  with  a volume  egual  to  that 

of  the  tr1ax1a!^ell1psoid  of  star  A (equivalent 

ratio  of  equivalent  radii  (a^jj/a^) 
gravity  exponent  for  ster  A 
gravity  exponent  for  star  8 
mass  ratio  (star  8 to  star  A) 
polytropic  Index  of  star  A 


Syrtol  Meaning 


polytrcplc  1nde»  of  star  B 
effective  wavelength  of  observations 
tcgarltftn  of  surface  gravity  of  star  A 
logarithm  of  surface  gravity  of  star  3 


For  systems  like  RS  CVn  In  which  the  stars  are  sensibly  spheri- 
cal the  egulvalent  radii  are  useful  parameters.  On  the  other  hand,  they 
have  very  little  physical  meaning  and  the  elllpsolOal  axes  are  to  be 

eters  (whether  It's  a spherical  radius,  a set  of  radii  for  a trlaxlal 
ellipsoid,  or  3 radii  of  a Roche  lobe]  Is  merely  a way  of  conveying 
quantitatively  information  about  a natural  phenomena,  which  most  likely 
does  net  actually  conform  (In  detail)  to  any  of  the  model  parameters. 

by  the  method  of  differential  least  sqaures  and  produced  a synthetic 
light  curve  [the  CALCM  dataset  referred  to  earlier).  The  integration 
was  performed  using  the  4xd  Gaussian  quadrature  option,  which  has  an 
accuracy  of  0.6G.  The  model  atmosphere  used  is  that  which  was  supplied 
by  Wood  In  Status  Report  #7. 

The  first  19  parameters  (ohly  16  of  which  are  used  here)  may, 
hy  user  selection,  be  improved  by  the  differential  least  squares  sub- 
routine. These  parameters,  for  the  purpose  of  this  research  were  div- 


Ug.  The  third-order  parameters  are  T^  and  QH.  The  fourth-order  param- 


The  other  parameters 


of  the  KINK  model  must  be  assigned  by  the  user  and  cannot  be  Improved 
by  the  differential  least  squares  routine. 

In  this  dissertation  the  terms  "free-parainetBP"  and  "adjustable- 
oarameter"  will  refer  to  a parameter  from  the  above  set  which  has  been 
designated  by  the  user  as  a parameter  to  be  Improved  by  the  differential 
least  squares  routine.  Likewise,  the  term  "fixed-parameter"  refers  to 
a paraamter  which  has  been  assigned  by  the  user  a particular  value  and 
not  allowed  to  be  changed  by  the  klINK  routines. 

No  third  light  was  allowed,  and  only  the  linear  limb  darkening 
law  was  used.  Certain  parameters  were  by  necessity  assisned  at  the  out- 
set and  never  altered.  These  parameters  and  their  assumed  values  are 


ASSUMED  PARAMETER  VALUES 


Parameter 

Assumed  Value 

j,  - I.B.,  a circular  orbit 

T^ 

67Q0"K  from  Popper  (1961) 

1 

1.D45  fron  Popper  (1961) 

log  g* 

3-5 

4.0  from  Gray  (19761 

log  ss 

3.5  fron  Gray  (1976) 

The  program  was  modified  to  allow  for  up  to  EGO  data  to  De  read 
or  written,  and  to  oroduce  an  additional  page  of  output,  which  is  merely 
a sunsnary  of  the  least  square  iterations.  In  all  other  respects  the 
program  Is  complete  through  Status  Report  fS,  The  version  of  tne  MINK 
model  used  here  was  updated  through  Status  Report  fS. 


The  ROHC  Prooram 

This  program  Is  very  simple.  It  reeds  the  CAL01  dataset  pro- 
duced by  UlIiK  and  converts  the  magnitudes  (m)  to  intensities  (I)  by  the 
relation 


This  conversion  is  performed  for  each  successive  record  in  the 

The  phase  of  the  two  records  are  checked  to  assure  that  they  are  the 
same.  If  they  are  not  the  same, an  error  message  is  written  and  the 
program  continues.  After  the  subtraction  of  the  CALC  record  intensity 
from  the  DBS  record  intensity  the  result  is  written  to  the  OHC  dataset 
along  with  the  phase  from  the  CBS  record.  These  calculations  are  re- 
peated until  all  the  records  have  been  processed.  It  should  be  noted 
that  the  CALC  dataset  is  written  for  plotting  by  the  LCPLOT  program 
later  in  the  procedure. 

The  FOURFIT  Prooram 

This  program  reads  all  the  records  from  CMC  dataset  and  by  the 
method  of  least  squares  determines  the  coefficients  of  the  following 
truncated  Fourier  series  which  best  represents  the  data  In  OMC. 

Since  the  data  in  and  around  secondary  eclipse  are  confused  by 
the  transit  of  the  smaller  star  across  the  assumed  non-uniform  surface 
of  the  larger  star  these  points  (from  phase  0.43  to  phase  0.57)  are 
excluded  frcm  the  least  squares  process. 

After  the  coefficients  of  the  truncated  series  nave  been  deter- 
mined they  are  used  to  evaluate  the  expression  at  each  value  of  phase  in 


trie  QSS  dataset.  The  result  of  each  evaTuatlon  is  trien  subtracted  from 
the  CBS  intensity  at  triat  phase.  The  results  are  what  one  would  expect 
If  RS  CVn  was  not  complicated  by  the  non-unifono  surface  brightness  of 
the  larger  star.  This  is  not  true,  hcuever,  for  the  secondary  eclipse 
because  the  photometric  effects  of  the  transit  geometry  cannot  be  so 
simply  removed.  These  effects  ^ contain  important  information  about 
the  surface  brightness  distribution  on  the  larger  star. 

The  other  functions  of  FOURFIT  are  to  produce  two  datasets  for 
plotting  by  LCPLOT.  The  first  dataset,  called  fOUC,  is  the  evaluation 
of  the  trisicated  Fourier  series  at  240  phase  points  sp  that  the  plot  of 
this  dataset  will  appear  to  be  a dontinuous  Tine.  The  second  dataset, 
called  TOMC,  is  the  difference  between  the  intensity  of  a record  in  the 
ore  dataset  and  the  evaTuatlon  of  the  Fourier  series  at  the  phase  for 
that  record.  In  otherworda,  it  is  the  difference  between  the  "real" 
distortion  wave  and  the  Fourier  fit  to  it.  These  points,  when  plotted, 
should  be  randomly  displaced  (by  only  a small  ^unt)  about  a straight 
line  of  slope  zero  if  the  elements  of  the  eclipsing  system  are  nearly 
correct.  Here  again,  the  region  of  secondary  ecliose  is  an  exception. 

The  LCPLQT  Prooram 

This  program  is  very  complex  and  for  the  present  purpose  need 
not  be  discussed  in  detail.  It  plots  the  datasets  the  user  wishes  as 

stricted  to  use  in  tiie  WAVE  procedure. 

Detailed  HAVE  Procedure 

The  detailed  JCL  WAVE  in-stream  procedure  appears  in  Fig.  9.  It 
differs  from  the  above  only  in  the  details  of  the  procedure,  but  not  in 


//SrSPRINT  oc 
//srsiN  00  01 

//SVSUTl  00  OuNMflt=K»l« 

//srsiffz  DO  0SN=1SPA8MS,SPACE»(TRK, (2,1). ROSE) .DISP-(NOH,P*SS), 
//  OC8-(RECPH-PB.LREa-SO,aLKSIZE«S4O0,0SORG»PS).UKIT=SYSDA 
//SEN2  EXEC  PGH-IEBGEHER,COKI>"(«,LT) 

//SrSPRINT  00  OUWt 
//SYSIN  CD  0«W« 

//SVSOrl  DO  DDNAHE-DATA 

//SYSUTZ  DO  D5N.S8CeS,SPACE"(TRK,(Z.l)BLSE),DISP=(NEVi.PASS), 

//  0CB»(RECm.FB,LREa>8O,BLKSIZE"64O0,DSORG-PS).UNlT=SYSDA 
//GEN3  EXEC  PGH>IEBGENER,C0N0>(4,LT) 

//SYSPRINT  OD  omm 
//SYSIN  CD  OlilWf 
//SYSUTl  DO  OONAHE-PRED 

//SYSOTZ  DO  DSfl=SSPRE0,SPAC£*(TRK,(2,l).RLSE],DISP-{NEW.PASS), 
//  0CB»(RECFM=FB,LRECL=80.B1.RSIZE"64[J0.DS0RG=PS).UHIT=SYS0A 
//GEN4  EXEC  PGK*IEBGENER.CSN0'(4,LT) 

//SYSPfilKT  DD  DUIW 
//SYSIN  DO  DUMMY 
//SYSUTl  DO  DDHAME'.  .. 

//SY5UTZ  00  DSK=«ST0P,SP4CE»(TRX.(Z,1),RI 
//  DCB=(R£CFH"FB,LRECL»80,eLKSlZE'6400.!En. 

//MNS  EXEC  PSM-IEBSENER,C0ND*(4.1T) 

//SYSPRINT  DD  DUWY 
//SYSIN  DO  Dim 

//SYSUn  OD  DSN'SSOBS.OISP*(OLD,PASS) 

OD  OSN'M5TOP,O!SP-(0LD,PAS5) 

//SYSUT2  OD  OSN-SSNOIW,SPACE"(TRK,(2,l),fiL! 

//  UNIT=SYSOA,DC8=BU!SIZE=8400,LRECL=M,RI 
//NINC  EXEC  PSH>NINX,C0ND’(4,LT) 

//FT05F001  DO  DSK-HPARMS,DISP-(OLO,PASS) 

//  DO  OSN-JWNKIN.DISP-iOLO.PASSl 

//  DO  DSH°8SPR£D,0ISP°10LD,PASS] 

//  DO  DSK-^tSNORM.OISP^OLD.PASS) 

//FTO6FD0I  DD  SVSOUT-A 

//FTD7FOO!  DO  OSK-&HINKOUT,SPACE-(TRK,  (2,1 ) 

//  DCB*(RECFH*FB.LRECL>8a.8LKSIZE'8400.0S( 

//FT08FD0I  00  SYS01FT"A 
//OMCA  EXEC  Pa=»MI*:,CQKD»(4,LT] 

//FTDIFOOI  DD  OSN=JSOBS,DISP=(OLD,PASS) 


//FTO3FM1  DD  OSN-SSDMC.SPACE'TRK. (2,1)  .RISE)  .0ISP"(N£W,PASS1 , 
n 0CS-(RECFM>FB,LR2CL>8O,BLKSlZ2>6«)O,0S0RI>PS).UNIT'SYSM 
//FTO6F001  OD  SYS0UT=A 


//  (e.d.rlse; 
Ntm  EXEC  PI 
//FTOlFOOl  or  • 


WRFIT,C0N(M4,LT) 


..  DC8>(RECFN>FB.LRECL*80.8lKSIZE-64O0,0S0ftS-PS),UIII>SVSM 
//PT06F001  DO  SYSOUT-A 

//FT18FQDI  DO  DSN-4&T0NC,0ISP-(NEV,PASS).l«IT-SYS0A. 

//  SPACE-(TM,{2.1).RLSE).DC8*(RECFM=FB,LREa=80,BLXSI2£*6400, 
//  DS0RS-P5) 

nm  EXEC  PS»-JEBSENER,CONO-{4,LT) 

//SYSPF 


//FH)8F001  DD  DSN*SJOeS,OISP-(OLO. DELETE) 
//FTO9F001  DD  DSN=SSC«.C,DISP»(OLD,DELETE) 
//FTlOFOOl  DD  DSrCitOMC.DlSP  [OLD. DELETE) 


//PTIDFOOI  DD  DSK=S&STDP.DISP  (OLD. DELETE) 
//PT14F001  DO  DSN»S4PARKS.0ISP  (OLD.OaETE) 
//PnSFOOl  DO  DSN"S4PR£0,0ISP  (OLD.DaETE) 
//FT16F001  DD  SCALCH.DSN  4SCALCM.DISP  (OLD. DELETE) 
//FT17F001  DD  DSN=S«R0eS.0I5P  (OID.DELETE) 
//FT18F001  DO  DSN=SSTOHC,DISP  (OLD. DELETE) 
//FT19FOOI  DD  DSN=S»FCHC.OISP  (OLD. DELETE) 

//PUNCK  EXEC  Pei*°lEBGENER.C0NO°<4,LT) 

//SY5PRINT  DD  DumY 
//SYSM  DO  DUMMY 

//SYSDTL  DO  DSM°tiROSS.DISP>(OLO, DELETE) 

//SYSUT2  00  SYS0irT-B,DC8-BUC5IZE*» 


Fig.  9.. 


ttie  HINK  parvneters  and 
model  atmosphere 


Che  085  dataset 

phase  and  intensity  read  hy  (2F1Q.5) 

/• 

//GEN3.PBED  DO  • 

the  "true-false”  card  (defines  Che  free-parameters} 

a|  o!o 

//GEK4.ST0P  00  • 

//LOT. PLOTS  00  * 

Che  LCPLOT  control  cards  for  the 
generation  of  plots  (see  fig.  8) 


Pig. 10. 'The  input  stream  which  uses  the  UAVE  procedure.  Note  that  a 
partitioned  load  module  library  must  be  supplied  with  the  necessary 


procedure 


86 

Sivei.  id  Fig. 

10  . The  missing  datasets  depand  on  the  application  of  the  procedune. 


CHAPTER  » 

AHALTSIS  AND  SOLLTriDN 
Introduction 

The  HAVE  procedure  eae  described  1n  the  preceding  diopter.  It 
W3S  develcped  to  assist  In  the  solution  of  the  light  curves  of  RS  CVn 
because  previous  techniques  were  Inadequate.  The  HAVE  procedure  Is  not 
capaole  of  arriving  it  i satisfactory  solution  without  the  critically 
iBiportant  supervision  of  a knowledgeable  human  being.  The  human  must  be 
capable  of  analyzing  Intermediate  results  and  determining  the  subsequent 
steps  by  wnidi  a better  solution  might  be  obtained.  At  this  oolnt,  HAVE 

not  yet  been  detennined  whether  the  procedure  can  be  applied  to  other 
eclipsing  or  non-eclipsing  systems.  It  1s  the  opinion  of  the  author  that 
HAVE  will  prove  to  be  helpful  In  obtaining  solutions, in  either  case, 
where  uneaplalned  photometric  complications  exist.  Therefore,  In  order 
to  assist  with  the  Implementation  of  HAVE  in  further  research,  a detail- 
ed account  of  the  process  by  which  the  present  solution  of  RS  CVn  was 
obtained  will  be  presented. 

The  Data  Analvzefl 

The  light  curves  of  RS  CVn  that  have  been  analyzed  with  the  use  of 
the  HAVE  procedure  are  listed  In  Table  8.  These  do  not  represent  all  of 
the  published  data  which  are  available,  but  they  are  a selectibn  which 
meets  the  requirements  of  the  research  effort  undertaken  as  the  topic  of 


Olssarcation. 


SOURCE  OF  DATA 


Publication  MOD  of  observations  Name 


Siturly  (19301 


IMS  Dissertation 
TJiis  Dissertation 
Tills  Dissertation 


5150A  1921 
SISOA  1963 

SISOA  1965 
SlSOA  1969 
SSOOA  197Sv 
aa90A  1975b 
3770A  1975U 
5SOOA  I97Bv 
4490t  1975b 


In  Table  3 the  last  coliaim  is  a "nacne"  which  has  been  given  to 
the  data  for  easy  reference  in  this  present  wort.  This  ‘‘name''  Is  also 
the  year  (and  color  in  some  cases)  in  which  the  data  were  obtained.  In 
a few  eases  It  will  be  convenient  to  group  the  1953,  1954,  1955,  and 


Catania  observatory. 

All  of  the  Catania  data  and  the  author's  data  of  RS  CVn  were 


obtained  using  BMS®  2420  as  a comparison  star.  This  greatly  facili- 
tates the  comparison  of  the  light  curves  frcoi  the  different  years, 
because  the  differentia)  magnitudes  are  in  the  same  light  units.  The 
difference  in  the  effective  wavelength  used  by  the  Catania  observers  and 
that  used  by  the  author  has  the  effect  of  changing  (only  slightly)  the 
effective  light  unit.  Since  the  bandwidth  of  the  filten  used  by  the 
Catania  observers  was  not  available.  It  was  imoossible  to  make  allowances 


for  any  {Hfferences  1o  tfte  two  instruneolal  syateis-  Fortunately,  to 
seme  extent  thfs  difference  does  not  cause  a oroblem  with  tne  analysis 
Secause  tlie  WINX  nde1  uses  the  effective  wavelength  as  an  inout  param- 
eter. Thus,  the  majority  of  the  HINK  paruaeters  that  are  subsequently 
determined  are  Independent  of  the  wavelength. 

In  Table  9 the  specific  dates  for  the  observations  are  given 
because  only  portions  of  the  published  data  were  used  for  the  1921  and 
Catania  light  curves.  The  migration  of  the  distortion  wave  would 
"wasn-ouf  some  of  the  detail  In  Che  light  curve  if  the  duration  of  the 
Observing  season  was  too  long.  For  this  reason  the  data  which  produced 
Che  light  curves  used  in  the  ilAVC  procedure  were  limited  to  those  points 
taken  over  a period  of  less  Chan  149  days.  The  data  taken  by  the  author 
was  deliberately  limited  to  a short  time  span  (65  days  in  1976  and  96 

The  modifications  to  HINK  (see  Chapter  IV)  extended  the  capacity 
Of  the  model  from  100  points  to  200  points.  This  would  have  allowed 
light  curves  to  be  made  with  up  to  200  normal  points.  However,  suoh 
light  curves  would  have  been  prohibitively  expensive  to  run  on  the 
computing  system  available.  Therefore,  a compromise  was  made  between 
the  cost  of  computing  (idilch  increases  approximately  as  the  square  of 
the  manber  of  normal  points)  and  the  desirability  of  a large  number  of 
nomals  which  would  give  a good  resolution  of  the  details  in  Che  light 
curves.  The  cempromised  value  was  150  normal  points  per  light  curve, 
or  as  close  to  this  value  as  was  realistic. 

The  calculation  of  the  normal  points  was  carried  out  on  the  dig- 
ital computep.  In  computing  the  normal  points  the  first  step  was  to 
cemputa  the  phase  for  each  individual  differential  magnitude  from  the 


linear  epfiemeris 

WO  {minlnnjBi  Hgtit)  " 4182S.042  t 4?797855E. 

The  DeriDd  of  this  eohemeris  was  coeiouted  Oy  Oacalano  and  Eodano  (1974). 
The  initial  e;iodh  Is  from  an  observed  lelnimun,  published  In  the  sene 
paper,  closest  to  the  start  of  Che  author’s  oOservatlons  in  1975.  Only 
the  declnal  fraction  was  retained  from  the  computation  of  tJie  phase. 

The  second  step  involved  sorting.  Into  phase  order,  all  magnitude-vs- 
phase  points.  This  process  was  done  separately  for  each  dataset  as  they 
are  listed  1n  Table  8.  The  final  step  In  the  nomal  point  calculations 
was  the  nest  important.  The  average  value  of  Che  phase,  and  the  magni- 
tude of  the  average  value  of  Che  intensity  for  all  points  within  many 
snail  phase  intervals  were  computed.  These  calculations  were  defined  by 
three  parameters.  The  first  two  parameters  together  defined  a section 
of  the  light  curve.  This  was  called  the  phese  range.  The  third  param- 
eter defined  a ’bin  siae'  within  the  phase  range.  For  eaample, 
suppose  that  a phase  range  of  0^0  to  1?0  1s  defined  (i.e.,  tne  entire 
light  curve),  and  a bin  siae  of  is  also  defined.  The  first  normal 
point  win  Pe  computed  from  the  Individual  points  between  phase  0.0  end 
0,1.  The  second  normal  point  will  be  computed  from  the  individual 
points  between  0?1  and  Q?2.  This  process  is  continued  until  the  last 
normal  point  is  computed  from  the  individual  points  between  0?9  and  1?0. 
In  this  enample  it  is  possible  to  have,  at  the  most,  10  normal  points. 
The  actual  number  would  depend  on  the  distribution  of  iridividuai  points 
throughout  the  light  curve. 

in  practice  it  is  better  to  have  more  than  one  phase  range  and 
bin  siae  tar  a light  curve.  There  are  two  reasons  for  this.  First,  the 
distribution  of  individual  points  will,  in  a11  probability,  be  unequal. 


Ttierefore,  it  would  be  advisable  to  cboose  the  phase  range  and  bin  aiae 
to  accomodate  the  real  distribution  df  data  in  the  light  curve.  Second, 
and  nore  importantly,  the  slope  of  the  light  curve  is  much  greater  during 
eclipses  than  in  the  rest  of  the  lipit  curve.  For  this  reason  it  is 
necessary  to  have  a snail  bin  size  during  an  eclipse  in  order  to  accurete- 
ly  represent  the  individual  points  with  the  nonral  points.  When  the 

considered.  The  final  set  of  normal  points  will  inevitably  be  a compro- 
siise.  The  data  which  are  ultimately  used  should  be  the  ones  which  best 
represent  the  individual  points.  The  nonnal  points  which  were  used  in 
the  WAVE  procedure  are  here  presented  in  Tables  12  through  22-  Notice 
that  each  table  gives  the  phase,  differential  mageitude  (variable  minus 
comparison),  and  nunCer  of  individual  points  used  for  each  hoimal  point. 

In  addition,  the  phase  range  and  bin  size  used  in  the  computations  ere 


There  were  many  computer  runs  using  the  WAVE  procedure  before  a 
final  solution  was  selected;  103  in  fact.  It  would  not  be  profllable 
to  discuss  each  of  these  runs  individually.  The  better  approach  is  to 
group  them  according  to  steps  in  which  the  runs  are  related  in  spme 
procedural  way-  This  reduces  the  nunber  to  approiimately  a dozen  steps, 
which  serve  to  slmpli^  the  discussion  which  follows.  In  the  future 
applications  of  the  WAVE  procedure  it  will  be  possible  to  reduce  the 
punter  of  runs  and  the  number  of  steps  considerably  based  on  the  experi* 


outline  of  the  process  by  which  the  light  curves  of  RS  CVn  were  solved 
vrill  be  presented. 

S swrrtery  of  the  solution  process  is  given  in  Table  9-  The 
step  names  are  listed  as  they  are  used  in  the  disoussion  which  follows 
The  values  for  the  parameters  in  Table  9 are  averages  of  the  computer 
runs  for  each  step.  The  standard  error,  below  each  value,  is  actually 
the  average  of  the  standard  errors  from  the  UlflK  program.  The  reader 


is  cautioned  that  tht 


error  is  listed,  • 


ire  computed  by  assuming  that  all  fixed-param- 
Thus  they  represent  an  optimistic  estimate  of 
the  parameter  was  determined.  If  no  standard 
ler  in  the  table  was  used  as  a fixed-parameter 


The  values  of  Table  i 
le  solution  process.  In  an: 
le  U light  curves  listed  ii 
eps  the  197Sv  ot 


are  only  a guide  as  to  the  evolution  o1 
given  step  it  Is  possible  that  not  all 
Table  8 were  run.  Furthermore,  in  manj 


light 


> light  curves,  or 
times.  The  1921  light  curve  is  c 
I averages  were  computed  by  giving 
;urves.  The  wavelength  o1 


precision  than  th 
different  fo 


averages  for  tl 


■es,  and  this  would  also  influence  the  validity  of  the 
.e  parameters  which  are  wavelength  dependent.  Therefore, 
lie  g are  to  be  considered  as  only  illustrative  of  the 
a final  solution  was  eventually  determined. 

At  the  outset  ef  the  solution  process  it  was  assumed  that  the 
solution  was  that  of  two  bodies  in  circular  orbits.  Therefore,  the 
eccentricity  of  the  orbit  and  the  value  of  third  light  were  assumed  to 


The  solution  process  wss  sUrted  0^  using  the  published  elensents 
[see  Table  2)  of  Popper  [19S1)  as  the  initial  paraeielers.  The  lS75v 
light  curve  was  used,  and  after  only  three  runs  a greatly  improved  dis- 
corcion  wave  resulted.  The  procedure  was  then  used  only  once  with  the 
1976V  light  curves.  These  four  runs  showed  that  the  WME  procedure 
would  give  improved  results  compared  with  previous  techniques.  This  step 
was  in  actuality  a preliminary  test  of  the  HAVE  procedure-  The  in-stream 
procedure  which  was  used  in  this  step  was  not  as  sophisticated  as  that 
which  was  developed  later.  The  in-stream  procedure  described  in  Chapter 


IV  was  used  in  steps  three  and  following. 

These  preliminary  results  showed  chat  Che  KO  star  corresponds,  in 
the  WINK  model,  to  star  8.  The  hotter  star  is  the  smeller  and  less  mas- 
sive star.  It  is  designated  by  the  'A'  subscript.  It  is  especially  Im- 


pertant  in  Che  case  of  R5  CVn  to  avoid 
and  "secondary”,  because  the  terms  are 


very  ambiguous.  For  example,  in 


RS  CVn  the  spectroscopic  primary  (i.e.,  the  more  massive 
the  photometric  secondary  ocmiponent  (i.e.,  the  fainter  star). 


meeting  of 
New  Nexico 

The  WAVE  procedure 


uoe  preliminary  in-sCream  procedure  was  used  in  step  two  as 
The  Catania  data  were  used  in  an  effort  to  obtain  plots 
low  the  evolution  of  the  distortion  wave  during  the  four 
by  these  data.  The  plots  were  presented  at  an  informal 
1 RS  CVn  working  group  of  lAU  Ccnifiisslon  42  at  Socorro, 
early  April  1978. 


i described  in  Chapter  IS 


and  first  used  in  step  three  to  improve  the  value  of  only  the  third- 
order  peremeters  (T^.  and  QH].  The  results  from  steps  one  and  tw  were 
preliminary,  but  of  sufficient  quality  to  be^in  a more  rigorous  process 
using  them  as  the  Initial  parameters.  The  third-order  parameters  are  not 
actually  elements  of  the  model,  but  they  must  be  determined  in  order  to 
arrive  at  a satisfactory  solution.  The  eleven  runs  of  step  three  were 
the  first  to  use  the  full  in-stream  WAVE  procedure.  Therefore,  this  step 
provided  the  needed  values  for  the  third-order  peranaters  and  the  SOBS 
datasets  necessary  for  subsequent  improvement  of  the  other  parameters. 
Since  T and  5M  are  unique  to  each  light  curve  no  meaningful  average  can 
be  computed  for  them.  For  this  reason  no  third-order  parameters  are 
Indicated  in  Table  9 for  step  three.  However,  it  is  of  interest  to 
ccmpare  the  average  HHS  error  for  the  11  light  curves  of  step  three  with 
the  subsequent  steps. 

Steps  Four  and  Five 

These  are  very  simller  steps.  They  both  represent  runs  which 
helped  the  author  gain  the  necessary  "feel"  for  the  effects  of  various 
parameters  and  cimnblnatians  of  parameters  on  the  solution  and  cost  of 
running  the  procedure.  Starting  with  step  five  the  value  of  log  gg  was 
changed  from  9.00  Co  3. SO  to  agree  with  the  published  values  ln  Oray 
(19:6). 

The  fourth-order  parameters  3/^  and  6g  were  made  free-parameters 
for  two  runs.  The  results  indicated  that  these  parameters  could  not  be 
Improved,  and  the  theoretical  von  Zeiple  values  were  adopted  for  both 
stars  (0.25).  Experimentation  with  the  value  of  o (mass  ratio)  and  the 
second-order  parameters  was  also  performed.  The  results  Indicated  that 
the  spectroscopic  mass  ratio  was  better  determined  Chan  the  photometric 


value  of  q,  arid  the  reflection  albedoes  were  very  poorly  determined  by 
the  data.  The  limb  darkening  coefficients  were  better  known. 


The  first-,  third-,  and  fourth-order  parameters  seemed  to  be  rea- 
sonably well  determined  based  on  prior  results.  Therefore,  It  was 
decided  that  the  second-order  parameters  should  be  Improved  before 
returning  to  Che  first-order  parameters  for  their  final  adjustment. 

This  was  the  expected  seouende  of  events.  However,  it  was  quick- 
ly found  that  the  solution  could  not  be  improved  by  having  only  the 
second-order  parameters  as  free-parameters.  It  was  necessary  to  add 
QM,  i,  and  T^  to  the  list  of  adjustable-parameters.  Various  combinations 
of  these  parameters  with  the  second-order  parameters  were  attempted  In 
order  to  Improve  Che  solution. 

At  the  conclusion  of  step  six  1t  was  felt  that  the  second-order 
parameters  were  as  well  determined  as  was  possible  with  the  available 
data,  and  they  were  never  changed  again.  The  values  which  were  adopted 
were  based  on  a selection  of  the  better  runs  of  all  previous  steps. 


This  step  consisted  of  only  three  ccmputer  runs  using  the  197dv 
date.  The  purpose  of  these  runs  was  to  establish  if  it  was  necessary  to 
improve  all  of  the  first-order  paremeters,  or  if  a smaller  subset  would 
be  sufficient.  It  was  hoped  that  the  inclination  and  temperature  of 
star  B were  sufficiently  well  determined  so  that  only  end  would 
need  to  be  imoroved.  Unfortunately,  it  was  found  that  all  four  param- 
eters had  to  be  changed  from  the  average  values  which  had  been 
In  order  for  the  HIHK  model  to  adequately  represent  the  ROSS  dal 


adopted 


From  tFio  results  of  step  seven  it  wes  obvious  that  the  average 
values  for  the  first-order  parameters  uere  not  adequate  to  represent  all 
of  the  light  curves.  So,  the  light  curves  for  all  the  datasets  in  TaOle 
8 were  run  with  all  four  first-order  parameters  as  adjustable-parameters. 
The  intent  was  to  find  a trend  In  one  or  more  of  the  parameters  as  a 
function  of  time  and/or  the  position  of  the  distortion  wave. 

No  clear  trends  appeared  evcept  for  an  unexpected  variation  from 
light  curve  to  light  curve  in  Che  value  of  the  Inclination.  It  seemed 
totally  impossible  that  the  Inclination  of  the  system  could  have  changed 
from  about  32^  In  1921  to  about  89^  in  1975.  This  Inexplicable  result 
led  td  the  sequence  of  runs  described  in  step  nine. 

The  1975V  data  were  about  the  best  of  any  used  in  this  analysis, 
so  this  light  curve  was  chosen  to  Investigate  the  nature  of  the  Inclina- 
tion of  the  system.  It  was  hoped  that  by  fixing  the  value  of  the  Incli- 
nation to  an  intermediate  value  (86?4]  between  the  extremes  found  in 
step  eight,  and  than  solving  for  various  combinations  of  the  ocher 
first-order  parameters  an  adequate  solution  could  be  found.  After  six 
computer  runs  It  was  deassed  impossible  to  find  a solution  with  an  incli- 
nation of  86?4  which  was  as  acceptable  as  the  one  previously  found  with 
1 ■ 89®. 


There  seemed  to  be  no  reesonable  explenatlon  to  th 
problem.  Therefore  a complete  review  of  previous  results 
see  if  some  error  had  been  made.  Indeed  an  error  had  been 


inclination 


s.  Originally 


a we11  known  value  because  the  color  of  this  star  could  be  observed  inde- 
pendent of  tfie  other  scar  during  primary  eclipse.  The  few  runs  of  step 
one  had  resulted  In  a temperature  of  SOSO^IC  for  star  A.  This  value  was 
then  used  as  a fixed-parameter  value  for  a11  subsepuenc  runs  and  T-  was 
adjusted.  While  this  would  not  cause  a large  error  In  the  detenclnaclon 
of  the  relative  brightnesses  of  Che  two  stars  It  was  noticeable  In  the 
determination  of  the  solution.  So,  step  10  was  a single  run  to  test  if 
taking  6700°K  (Popper  1961}  as  the  value  of  T,  could  resolve  the  Inclina- 
tion probiem. 

As  one  might  have  expected,  the  temperatures  of  the  stars  had 
little  affect  on  the  value  of  the  Inclination.  However,  it  did  Improve 
the  Rf^  error  of  the  WIWK  determination  for  the  1975v  light  curve. 

Ac  Che  conclusion  of  the  run  It  was  realized  Chat  Che  inclination 
problem  was  not  really  a problem  at  all.  It  was  merely  a msnifestatlon 
of  the  non-un1fonn  surface  brightness  of  the  KO  star.  Since  the  struc- 
ture of  this  surface  brightness  distribution  1s  a changing  function  of 
time,  the  transit  of  the  smaller  star  will  mimic  different  Inclinations. 

The  normal  points  In  secondary  eclipse  were  removed  from  all  Che 
light  curves.  Using  the  resulting  datasets  as  input  OBS  decks  a new  set 
of  runs  were  made  with  1,  T , t}M,  T_,  a .,  and  k as  free-paranetars. 

The  new  values  for  these  first-order  paranetars  were  remarkably  similar 
fpr  each  light  curve,  for  example,  the  range  in  the  value  of  the  Incll- 
natlpn  was  from  84?3  for  1964  to  86?1  for  1975u. 

The  results  from  step  II  almost  completed  the  solution.  There 
was,  however,  one  other  use  for  the  WAVE  procedure  which  could  add 


Insight 


pulsating  1 t would 


teaparatura 


and  radius  of  this  star  would  be  different  for  each  light  curve.  Chi  the 
other  hand,  1f  the  surface  of  the  KO  scar  is  hon-uji1fortn,  the  radius 
would  not  change*  only  the  ceinperature  would  vary. 

Kith  the  above  conditions  in  mind  and  the  points  of  secondary 
eclipse  still  removed,  another  set  of  11  runs  was  made  with  only  Tg  and 
h as  free-parameters.  The  radius  of  star  A was  a flMd-paraneter, 
so  the  changes  in  k would  be  solely  due  to  the  radius  of  star  B ta^gj- 

The  value  of  a^j  was  almost  identical  for  all  light  curves,  while 
the  value  of  was  somewhat  different.  This  would  support  a starspot 
mode!  over  a radial  pulsating  model. 


Sup  Thirteen 

The  last  step  was  to  return  the  secondary  eclipse  points  to  the 
datasets,  and  run  the  WAVE  procedure  Co  proOuce  a new  set  of  HOBS  data- 
sets. The  solution  portion  of  KINK  was  not  used  because  the  solution 
from  steps  11  and  12  was  felt  to  be  the  best  that  could  be  achieved. 
This  step  was  Intended  only  to  provide  the  final  ROBS  and  OfIC  datasets. 
Table  9 does  not  nave  a listing  for  this  step  because  there  were  no 


changes  from  step  12. 


il  Solution 


It  was  noted  In  step  12  above  that  T.  was  the  only  element  which 
s significantly  different  In  the  11  light  curves  which  have  been  solved 
tn  the  WAVE  procedure.  The  value  of  Tg  which  has  been  determined 
ir  each  light  curve  is  an  average  over  the  non-uniform  surface  of 
le  KD  ster.  In  addition  to  different  values  of  T.,  the  various  light 
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curires  also  have  differentFoorUr  coefficients  to  represent  tne  dis- 
tortion waves.  Table  10  ^Ives  the  elements  of  the  system  which  are 
cormon  to  a11  lioht  curves.  Table  II  gives  the  values  which  are  unique 
for  each  light  curve  as  determined  from  steps  11,12  and  12. 

The  errors  quoted  in  Table  10  are  standard  errors  determined  with 
the  use  of  Che  HAVE  procedure.  The  correlation  among  the  systmn  parame- 
ters makes  It  difficult  to  assess  the  precision  of  the  determination  of 
each  value.  The  Ideal  method  to  obtain  an  estimate  of  this  precision 
would  be  to  use  all  the  HIKk  model  parameters  as  free-parameters  and 
thus  allow  the  least  squares  routine  to  determine  the  standard  error. 

The  only  method  currently  available  with  the  HAVE  procedure  technique 
Is  to  use  Che  maximum  number  of  free-parameters  allowed  (9}  to  estimate 
the  stendard  errors.  By  repeating  this  two  additional  times.  It  Is 
possible  to  include  all  of  the  parameters  which  can  be  adjusted.  This 
was  done  using  the  lS75v  Tight  curve.  The  standard  error  given  in  Table 
10  is  the  maximim  value  of  the  standard  error  from  the  three  runs. 

The  Distortion  Have 

The  Use  step  of  the  solution  process  produced  a set  of  Fourier 
coefficients  which  represent  the  distortion  wave.  These  coefficients 
are  given  In  Table  11.  The  standard  error  of  each  coefficient  is  given 
below  the  parameter  value.  As  a reminder,  and  for  the  convenience  of  Che 
reader,  the  truncated  Fourier  series  from  which  these  coefficients  have 


The  Ap  term  is  a constant  which  resolves  the  nonnallzacloh  differences 
between  the  OBS  light  curve  and  the  CALC  light  curve.  Fresioiably,  had 
1t  Been  possible  to  accurately  normalUe  the  observed  light  curve  to 


unity  outside  eclipse,  QH  and  A would  both  be  zero.  Unfortunately, 
this  was  Impossible  to  do,  a priori.  In  the  case  of  RS  CVn  because  of  the 
large  asynnetrlc  and  dynamic  distortion  wave,  as  well  as  the  large  period 
changes  in  the  systen. 

The  other  terms  of  Che  truncated  Fourier  series  represent  the 

tenn,  and  the  apparent  nonnallzed  luminosity  of  star  B (Ig,  a quantity 
printed  by  WINK)  to  produce  a light  curve  of  the  KO  star  as  1t  would  ap- 
pear if  the  brighter  star  were  removed  from  the  system.  If  the  OhiC  data- 
set is  used  In  place  of  Che  truncated  Fourier  series,  the  result  will 
include  the  original  observational  error  which  has  been  partially 
Smoothed  from  Che  data  by  use  of  the  truncated  Fourier  series. 

The  computations  which  are  outlined  above  were  performed  on 
each  light  curve.  The  equation  that  represents  this  process  is 
(7)  Lj(e)  . iwc(e)  + 

where  Lg(6)  1s  Che  intensity  at  phase  e of  the  Isolated  K0  scar,  OMC  (a) 
is  the  Intensity  of  the  distortion  wave  at  phase  S , Ig  Is  the  nonnallzed 
apparent  luminosity  of  star  B,  and  QM  is  the  quadriture  magnitude.  The 
latter  two  quantities  have  been  listed  for  each  light  curve  in  Table  11. 
The  above  equation  produced  light  turves  of  the  KO  star  which  are  ln 
units  of  the  comparison  star's  brightness.  These  light  curves  are  dis- 
played in  Figs.  11  through  21.  The  region  of  secondary  eclipse,  in 
these  figures.  Is  conpl icatod  by  the  transit  phencmenon,  and  therefore 
this  region  does  not  represent  the  light  changes  due  solely  to  the  KO 

The  values  of  distortion  wave  minimum  and  maximum  and  6,ggj,J 
ape  listed  in  Table  11.  They  have  been  computed  by  differentiating  the 
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Fourier  series  of  equation  {6}  witti  coefficients  frcn  Table  11,  setting 
tne  result  equal  to  aero,  and  evaluating  the  solution  by  Neivton's  oettiod. 

Period  Changes 

If  the  period  of  a binary  system  changes  vrith  time  it  is  an 
Indication  that  the  angular  nomentian  of  the  systen  has  changed  or  that 
the  distance  to  the  system  has  changed.  There  are  many  ways  In  i^1ch 
Che system^sangular  momentum  can  be  altered.  The  one  which  con«s  to  mind 
most  frequently  is  the  transfer  of  mass  from  one  star  td  the  dther 
through  the  inner  lagrangian  point.  This  method  of  mass  transfer  has 
been  of  great  Interest  during  the  lest  two  decades  as  a means  of  eaplain- 
tng  the  evolution  of  binary  stars.  It  is  mentioned  here,  not  as  a 
possible  explanation  for  period  changes  In  RS  CVn,  but  as  a reminder 
that  it  is  only  one  of  several  ways  to  change  the  angular  ibcmentum  of 
a binary  system,  and  only  currently  enjoys  a large  degree  of  attention. 

This  is  not  intended  to  miniraiae  the  importance  to  astronomical 
knowledge  of  any  particular  period  change  mechanism,  but  rather  to  help 
explain  the  current  state  of  affairs.  Ratten  (1973)  has  a good  review 
of  period  change  mechanisms,  briefly  they  are  as  fallows;  apsidal 
motion,  the  presence  of  e third  body,  and  real  changes  in  orbital 
period  which  maybe  due  to  mass  transfer,  mass  loss,  or  mass  gain. 

The  later  three  mechanisms  are  in  actuality  changes  in  the  angular 
momentum  of  the  syston.  The  mass  exchange  due  to  evolution  Is  simoly 
a special  case.  Of  course,  there  is  always  the  possibility  of  another 
unknown  mechenlem. 

R3  CVn  and  RS  CVn  binaries  may  prove  to  be  extremely  important 
systems,  because  it  Is  reasonably  certain  that  the  period  chenges  are 
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not  doe  to  a 11gnt-t1me  effect  (orbital  mollon  aOoot  a third  body)  or 
to  mass  tranafer-  Keither  Is  It  lihely  that  mass  loss  nor  apsidal 
motion  are  the  causes  of  the  ouasi-periodic  period  changes.  The  research 
in  the  coning  decades  nay  return  to  one  of  these  causes  of  period 
change,  but  at  the  present  all  seen  Impossible.  Therefore  it  Is 
important  that  the  period  be  measured  as  accurately  as  possible.  In  the 
case  of  star  systems  like  RS  CVn  which  have  very  asyimieCric  primary 
eclipses  lt  Is  difficult  to  accurately  determine  the  time  of  minimum. 

This  15  compounded  by  the  long  periods  (and  earth-bound  photometry). 

The  iilHK  parameter  T^  can  be  determined  to  an  accuracy  equivalent  to 
less  than  40  seconds  in  the  case  of  RS  CVh.  Unfortunately,  tnis  is 
somewhat  misleading  because  the  absolute  time  is  not  knowm.  The  para- 
meter T is  that  which  best  represents  a light  curve.  The  mean  epoch 
for  the  data  which  comprise  the  light  curve  Is  difficult  to  determine 
because  the  computation  is  complicated  by  the  random  nature  of  the  data 
acquisition  and  Che  unequal  effect  qn  T of  eaqh  datum.  Nevertheless, 
the  parameter  can  be  a useful  tool  in  deCennlning  the  nature  of 
period  changes  since  It  affords  a higher  degree  of  confidence  than 
other  methods.  This  trill  allow  the  investigation  of  angular  imnehtup 
changes  due  to  processes  that  are  more  insidious. 
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CHAPTER  VI 
DISCUSSION 


Intro  ductiori 

The  topic  of  this,  the  lest  chapter,  will  be  the  impi icetlone  of 
the  data  presented  in  tiie  foregoing  five  chapters.  First,  the  properties 
of  the  distortion  wave  will  be  considered.  The  dates  of  TiBxiinueT  and  orin> 
inusi  Tight  will  be  competed.  The  resulting  goantities  will  be  referred  to 
as  tnd  respectively.  From  these  values  it  wil  1 be  possible 

to  coicpare  the  peab-to-peax  amplitude  for  each  of  the  II  Tight  curves 
under  investigation.  The  vaTues  for  the  phase  at  minimum  and  the  ampTi- 
tude  of  the  distortion  wave  can  then  be  compared  with  those  in  Kali  U972, 
197Sa).  Sy  taking  advantage  of  the  compTication  in  the  distortion  wave 
around  secondary  eclipse  due  to  the  non-uniform  intensity  of  the  surface 
on  the  KO  (star  B)  component  an  estimate  of  the  teimierature  of  the  star- 
spot will  be  made.  The  discussion  of  the  distortion  wave  wiT1  conclude 


with  some  doninents  about  the  nature  of  thi 
TTie  problems  associated  with  the  r 
be  considered.  The  existence  of  short-ter 
pointed  out  by  KalT  (1972,  1975a),  are  not 
and  realistic  KAVE  procedure  model.  The  n 
remain  unexplained,  at  least  in  classical 


on-uniform  period  changes  will 
n period  variations,  first 
confirmed  by  the  more  accurati 
uch  longer  period  variations 


The  Distortion  Wave 

The  value  of  the  phase  at  meximum  and  itrinimimi  Tight  of  the 
distortion  wave  are  listed  in  Table  11.  It  would  be  more  convenient  if 
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tti«se  val  Lies  i^re  not  dependent  en  tne  epneineris  uhlcti  was  used  Co 
udinpuce  tlie  ^eses  fon  the  construction  of  the  11gnt  curves.  As  a first 
step  in  this  direction,  nev  estimates  for  Che  tires  of  minimum  (hJD^) 
were  compuCed  for  each  light  curve.  These  values  for  HJD  are  listed 
in  Table  23.  They  were  computed  by  using  Che  value  of  T (from  Table  11) 
Co  correct  the  time  predicted  by  the  linear  e^emeris  ^equation  4). 

The  published  values  of  the  times  of  minimum  (Catalano  and  Rodono  1974) 
for  the  Catania  data  were  used  as  a guide  In  establishing  the  number  of 
oroital  cycles,  for  the  light  curves  of  1975  and  1975  a value.  ft)6, 
was  computed  by  averaging  the  observed  dates  of  primary  mlnimian.  weighted 
by  the  square  root  of  the  nunber  of  points  on  each  date.  This  value 
was  then  used  to  coreute  a cycle  nianber.  and  thus  a value  of  liOD  for 
each  of  these  six  light  curves,  by  Incorporating  tbe  appropriate  values 


The  next  step  1n  computing  estimates  for  the  dates  of  maximum 
and  minimum  light  of  the  distortion  wave  and  hldO^j,)  Involved 

the  calculation  of  the  mmber  of  days  after  MJO^  chat  these  events 
occurred,  from  Table  11  the  values  of  and  can  be  obtained. 

By  multiplying  these  values  by  the  period  (4?7976S6)  used  in  the  linear 
ephemeris.  the  appropriate  delay  (In  days)  was  obtained.  This  would 
have  been  somewhat  inaccurate  because  and  from  Table  11  are 
dependent  on  tne  linear  ephemeris.  To  reduce  this  effect,  the  appropri- 


values  that  were  computed  by  th 
ly  published  values. 


10  listed  in  Table  23.  The 
sore  reliable  than  previous- 
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Rg.  22. -The  visual  band  aiqilltude  of  We  dlstor- 
tion  wave  olotted  as  a function  of  the  ortiltal 
cycle  (E)  from  Hall  [19721. 


Hall  predicted  Wat  We  variation  in  auplitude  of  We  distortion 
nave  xas  We  result  of  a 23Vyear  "sotisoot  cycle"  operating  in  We  sub- 
giant  of  the  system.  Fig.  22  shows  a piot  of  distortion  wave  amplitude 
in  We  visual  band.  The  filled  circles  are  from  Hall's  (1972)  paper, 
while  We  crosses  are  from  Table  23.  The  abscissa  is  plotted  using 
Hall's  emphemeris.  The  solid  line  is  We  functional  dependence  pre- 
dicted by  We  1600  orbital  cycle  [23Vyear)  "sunspot  ^cle."  Note  Wat 
We  new  observations  are  In  direct  conflict  with  We  predictions,  but 
We  new  estimates  of  We  amplitude  from  the  Catania  data  agree  very  well 
wlW  Hall's  graphical  astimates. 

This  1s  ciroumstantial  evidence  Wat  We  phenomenon  In  HS  CVn  is 
not  We  same,  in  all  respects,  as  Wat  known  as  sunspots.  2oW  effects 
are  probably  surface  phenoinena  resulting  in  a non-unifora  surface  bright- 
ness. but  Wey  are  quite  different  in  oWer  ways. 

AnbWer  consequence  of  Kali's  "sunspot"  model  was  a periodic 
variation  in  We  migration  rate  of  We  distortion  wave.  The  data  for 
1975  and  1976  are  again  in  conflict  wlW  We  evtrapclation  of  Wis 
relationship.  Hall  plotted  We  phase  of  distortion  wave  nrlnimum  as 
a function  of  orbital  cycles.  The  difference  between  We  observed 


E > 0 


more  ^et^i1.  Sfnce  the  analogy  with  the  sansoot  has  not  withstood 
observational  testing,  a more  detailed  discussion  will  not  be  given. 

The  distortion  wave  in  RS  CVn  is  not  the  result  of  a phencnenon  similar 
In  a11  respects  to  sunspots;  at  least  the  similarity  Is  not  as  great 
as  It  was  originally  portrayed. 

The  distortion  wave  Is  almost  certainly  oaused  by  temperature 
variations  on  the  surface  Of  the  cooler  star  of  the  system.  If  the 
size  of  the  cooler  region  Is  larger  than  the  disk  of  the  smaller  star 
as  the  transit  Is  viewed  from  earth,  then  an  estimate  for  the  temperature 
In  the  cooler  region  cen  be  mede.  In  addition,  the  amplitude  of  the 
distortion  wave  also  gives  an  estimate  of  the  temperature  variations. 

The  latter  approach  will  be  considered  first. 

The  average  temperature  of  star  6 (Tg]  and  the  nonnalized 
apparent  luminosity  (1g)  can  be  used  to  estimate  the  value  of  a in  the 
following  relation: 

(8)  1,.10-0-^‘i^warM. 


give  an  average  value  of  1.6  * 10’^^  for  a fin  units  of  the  comparison 
star's  energy  output  per  unit  orbital  radii  squared  per  degree  Kelvin  to 
the  fourth  power),  The  maaimum  temperature  on  the  KO  star  can  be  esti- 
mated by  assijnlng  that  the  hemisphere  facing  Che  earth  at  the  maximLin 
of  the  distortion  wave  is  not  spotted.  This  value  is  approximately 
SOOO'k,  with  a range  from  4900“k  (1963)  to  5100“s  (1975v).  Simllerly, 
the  temperature  of  the  cooler  region  can  be  estimated  to  be  about  4400^K. 

Another  way  to  estimate  the  maximum  tanperature  of  star  6 is  to 
assume  Chat  the  hemisphere  facing  the  earth  during  secondary  eclipse  In 


1975  ana  1976  »as  the  brighter  unspotted  hemisphere.  The  smaller 
star  will  then  black  more  light  then  the  UINK  model  predicted  for  star  S 
with  an  average  temperature  of  T^,  This  method  gives  a temperature 
of  approximately  SlOO^h.  This  Is  the  same  as  the  previous  method  and 
would  therefore  give  the  same  temperature  for  the  cooler  region  as 


The  temperature  of  the  spotted  region  can  be  estimated  by 
making  similar  asstanptions  as  were  made  above.  In  actuality  the 
estimates  are  only  upper  limits  on  the  temperature  of  Che  cooler  region. 
The  upper  limit  on  the  tesnperaCure  for  a spot  covering  one  hemisphere 
1s  Ifrwii  the  brightness  Of  star  5 at  distortion  wave  minimimi,  as  was 
ccmputed  above)  44C10^K.  The  probing  effects  of  Che  smaller  star  during 
transit  can  be  put  to  advantage  only  for  the  1963  light  curve.  In 
this  case  the  upper  limit  for  a spot  larger  than  the  apparent  disk  of 
Stan  A is  44O0”K.  The  similarity  in  value  from  the  two  estimates  would 
Indicate  that  the  size  of  the  active  regioh  is  almost  the  size  of  one 
hemisphere  of  star  B.  Ih  addition,  it  indicates  that  the  spotted  area 
is  probably  syitmetrldai  with  the  eQuacor  of  star  B,  unless  it  Is  so 
asymmetric  as  to  be  entirely  on  one  side  of  the  eguator  of  star  8. 

Before  leaving  the  subject  of  spots  it  is  important  that  the 
reader  differentiate  between  the  spots  of  RS  CVn  and  those  of  the  sun. 


it  conveys  the 


on  the  sun  may  be  quite  different  from  that 
pt'*  has  been  used  because  It  1s  convenient  and 


a hotter  photosphere.  The  shape  of  this  region  Is  not  circular,  neces- 
sarily, This  1s  demonstrated  by  the  shape  of  the  distortion  wave  In 
the  1963  light  curve.  The  shape  of  the  spot  on  the  KO  star  would 


be  SQo«thlTi3  like  thet  shown  in  Fig.  24.  This  shape  is  indicated 

sine  wave.  Since  the  distortion  wave  is  never  fiat  at  maximum  In  any 
of  the  light  curves  studied  so  far  it  is  necessary  that  the  sootted 
region  extend  at  least  ISO  degreesin  longitude  on  the  star. 

It  is  also  apparent  from  Che  double  miniimsn  in  1976  (see  Fig. 

17  in  Chapter  V]  that  the  concentration  of  spots  is  not  limited  to  the 
region  which  migrates  (approximtely)  uniformly.  The  significance  of 
this  is  not  clear,  but  certainly  any  model  devised  to  explain  the 
photometric  distortion  wave  must  be  capable  of  reproducing  the  rapid 
changes  in  the  detailed  shape  of  the  distortion  wave.  This  is  best 
visualiaed  by  comparing  the  changes  in  shape  from  197S  to  1976  (see 
Figs.  16  and  17),  or  the  changes  from  1963  and  1964  (Figs.  12  and  13). 

The  radial  pulsation  of  star  B cannot  be  rasposible  for  the 
photometric  distortion  wave.  If  the  pulsation  of  starBis  assumed,  Chen 
It  would  cause  a change  in  system  brightness  and  a subsaquent  change 


cause  ttie  observed  ampHtude  of  tbe  distortion  wave>  the  radios  would 
nave  to  change  by  approaltnately  20%.  The  data  studied  here  represents 
only  about  a quarter  of  a migration  cycle,  thus  a change  of  about  5S 
should  present  Itself  if  radial  pulsation  is  tne  cause  of  Che  distortion 
wave-  The  variation  in  a^j  in  tne  computer  runs  of  step  11  in  Chapter  V 
show  a 2S  variation,  at  the  most.  Additional,  and  higher  quality  data 
could  improve  the  confidence  with  which  radial  pulsation  is  rejected  as 
a cause  for  the  distortion  wave.  However,  the  results  here  are  reason- 
ably convincing  that  radial  pulsation,  from  observational  data,  is 
not  responsible  for  this  anomalous  photcmetric  behavior. 

Obscuration  of  the  light  from  the  IB  star  by  intervening  materiel 
as  an  axplanation  for  the  distortion  weve  has  been  rajected  by  Oliver 
(1974,  1979).  Other  possible  explanations  include  non-radial  pulsation, 
unusual  convection  patterns,  magnetic  effects  (as  in  sunspots),  and 
"soraething  else."  Time  will  perhaps  eventually  show  what  is  the  correct 
axp1 anation. 

Qrbital  Period  Variations 

Tha  short-tarm  period  verietions  in  RS  CVn  have  been  attributed 
to  mass  ejedtion  by  Hall  (1972,  l97Sa).  Hall  claimed  that  mass  was  being 
ejected  frcm  the  brighter  hemisphere  of  star  B.  When  this  hemisphere 
was  feeing  in  the  direction  of  orbital  motion  it  would  tend  to  reduce 
Che  period.  In  Fig,  25  are  plotted  the  times  of  minima  as  corrected  Py 
Hall  (1975a)  for  the  effects  of  the  distortion  wave,  and  the  times  of 
minima  as  they  appear  in  Table  23.  Kail's  estimates  ere  the  filled 
circles,  while  those  from  Table  23  appear  as  crosses. 

The  data  in  Fig.  25  are  from  Catalano  and  Rodono  (1974),  and 
wera  published  as  times  of  minima.  Hell  corrected  tnese  values  by 


electric  tlines  of  mlniraa-  Pilled  circles  ere 
points  corrected  for  the  distortion  «ve  effects 
Py  Hall  (197Se).  Crosses  are  points  as  cotoputed 
in  this  dissertation.  Arrows  Indicate  when  the 
period  is  predicted  to  decrease  and  Increase 
according  to  Hall  (1972,  1975al. 


estimating  the  effect  of  the  distortion  wave  on  the  actual  times  of 
conjunction.  The  HAVE  procedure  has  accomplished  the  same  thing  (albeit 
with  more  effort)  with  greater  accuracy.  Unfortunately,  it  was  only 
possible  to  do  this  with  the  published  data  for  1963,  196A,  1965,  and 
1966,  and  of  course.  Che  1975  and  1976  data  obtained  by  the  present 
writer.  The  arrows  Indicate  idien  Hall  predicts  an  increase  (or  decrease) 
In  orbital  period  due  to  the  mass  ejection  from  the  brighter  hemi- 
sphere. note  that  with  the  change  of  the  single  point  for  1966  the 
variation  bectmies  less  obvious.  In  addition,  the  predicted  decrease 
aoparently  did  not  occur  in  1975  or  1976.  It  would  appear  that  Hall 's 
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It  is  difficult  to  iTsgine  ttat  the  short-term  variations  whlcn 
apcear  to  renain  are  observational  error.  This  is.  ho«ver,  a possi- 
bility. If  they  are  Indeed  real,  then  an  adequate  esplanation  would 
probably  follow  along  similar  lines  to  those  of  Hall,  tut  with  less  reli- 
ance on  the  sunspot  analogy.  On  the  other  hand,  it  would  be  expected 
that  both  stars  are  losing  mass  r1a  a stellar  wind,  since  they  are  very 
similar  1n  age  and  mass.  If  this  is  the  case,  then  variations  In  the 
difference  In  the  loss  rates  might  cause  the  short-term  period  changes. 
The  data  presented  here  do  not  show  that  any  short-term  variations  exist. 
If  they  do  not,  the  mechanism  for  the  much  longer  term  variations  also 
vanishes.  Here  again,  a stellar  wind  nay  Pe  operating  In  one  or  both 
stars  to  cause  tbe  (aimostl  periodic  0-C  variations.  If  this  cannot 
account  for  the  variations,  then  It  may  be  necessary  to  return  to  the 
consideration  of  a very  under-luminous  third  body,  such  as  a white  dwarf, 
neutron  star,  or  black  hole.  The  latter  two  possibilities  are  highly 
unlikely, 

Cone! usioh 

Additional  light  curves  of  RS  CVn  have  been  bbtained  as  a result 
of  this  research.  They  are  the  first  complete  light  curves  ih  more  than 
one  color  to  be  published.  The  light  curve  for  1976  hes  shown  a pre- 
viously unknown  characteristic  of  the  system,  namely,  that  the  distor- 
tion wave  can  consist  of  more  than  one  mlnlimim.  The  1976  distortion  wave 
also  shows  an  amplitude  as  large  or  larger  than  any  previously  published. 

The  computer  progrem  UINK  and  associated  programs  which  coerprised 
tbe  HAVE  procedure  have  previded  a sec  of  elements  which  are  better  than 
any  previously  obtained.  The  technique  has  also  provided  a mathematical 
characterization  of  the  distortion  wave  which  can  be  treated  by  less 
complicated  models  then  are  currently  being  used  (Eaton  1977).  This 


s’tmpMflc&t'lDn  1s  due  Cd  trie  less  coiaplicated  case  of  studying  t^e 
phenomenon  of  the  distortion  wave  with  the  "complication''  of  the  eclipses 
remved,  or  nearly  reenved- 

The  additional  light  curves  and  improved  solution  method  have 
added  moTe  weight  for  the  rejection  of  radial  pulsation  as  the  source  of 
the  distortion  wave.  Popper  (1961)  had  far  too  little  information  with 
which  to  investigate  Che  question.  CataTanc  and  Rodono  (1966)  perhaps 
had  sufficient  data,  but  ttiey  lacked  a better  and  more  efficient  method 
of  solution.  Hall  (1972)  was  in  a similar  predicament.  Now,  however,  e 
measurement  of  the  radius  of  star  Q to  sufficient  accuracy  1s  possible. 
The  analysis  of  the  data  has  shown  that  the  needed  verlatlon  In  r^  Is 
not  seen.  Confirmation  of  this  with  additional  data  1s  desirable  In 
order  to  confirmi  that  radial  pulsation  Is  not  responsible  for  the  distor> 

The  constancy  of  Che  comparison  star  61>r35^  2420  has  been  con- 
firmed (sea  Fig.  6). 

The  analogy  between  Che  nonuniform  surface  brightness  on  the  KO 
star  of  RS  CVn  and  classical  sunspots  has  been  shown  to  be  over-rated. 

The  amoTItude  of  the  distortion  wave  does  not  appear  to  follow  a quasl- 
penodlc  vanatlon  as  closaly  as  does  the  nirmber  of  sunspots  on  the  sun. 
Furtnemore,  the  nearly-perlodlc  variation  1n  the  migration  rata  of  the 
distortion  wave  did  not  continue  into  1975  end  1976  as  was  predicted. 

The  nngraClon  rate  was  nearly  constant  from  1963  through  1976,  showing  a 
slight  decrease  1n  rate  rather  than  the  predicted  increase.  Finally, 
the  short-term  period  variations  which  were  considered  an  integral  part 
of  Che  sunspot  analogy  were  not  seen  in  the  data  used  in  this  Investiga- 
tion. In  fact,  a new  ephemerls  can  Pe  computed  for  Che  Catania  and  1976- 
76  data  studied  here.  The  resulting  ephemeris, 


ISl 

(9)  WJO,  “ 42870.962  ♦ 4?7977S2E. 
indiutes  that  the  orbital  period  fron  1963  to  1976  leas  pearly  coPStant. 
The  naxlmim  deviation  from  the  observed  times  of  minina  in  Table  23  is 
6 minutas  CO. 004  days)  for  the  1966  light  curve. 

The  shape,  temperature,  and  life-time  of  the  spots  have  beep 
estimated.  The  shape  of  the  distortion  wave  is  frequently  a "saw  tooth" 
which  implies  a triangular  or  tapered  shape  Csee  Fig.  24).  the  upper 
limit  for  the  temperature  of  the  spotted  region  has  been  set  at  44Q0®K. 

The  development  time  for  a large-scale  spot  must  be  approximately  one 
year.  The  smaller  minimum  in  the  distortion  wave  was  clearly  visible  in 
1976,  but  there  was  only  a hint  of  its  existence  in  1979.  In  addition, 
an  estimate  of  the  surface  temperature  on  the  unspotted  photosphere  was 
made;  it  was  9100*K. 

Additional  research  is  needed  before  a complete  model  for  RS  CVn 
can  be  developed.  High  time-resolution  photometry  and  apectroscopy  in 
the  optical  region  is  desirable  and  would  be  especially  useful  during 
eclipses.  Secondary  eclipse  could  yield  more  information  about  the  nature 
of  the  spots,  i.e.  unifonn  versus  blotchy  spots  and  a more  detailed  lup 
of  the  distribution  of  the  spot.  Primary  eclipse  could  provide  informa- 
tion on  the  atmosphere  of  the  M sUr  from  observations  during  the  four 
contact  points. 

Of  course,  observations  in  other  spectral  regions  would  provide 
needed  infonnation.  In  addition,  more  information  from  the  X-ray  and 
radio  regions  would  be  helpful  in  determining  If  the  distortion  wive 
phenCBtienon  is  associated  with  the  emission  at  tbesa  wavelengths. 


binaries. 
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ei0eun<!ui  SKETCH 


On  17  Oecenter  1949,  Elwyn  Whit  Ludington  was  bom  to 
Sergei  A.  and  Laona  Inez  Ludington  in  Eugene,  Oregon.  He  attended 
public  school  in  Eugene  and  spent  the  suomer  months  of  his  teenage 
years  with  his  uncle  and  aunt  (Robert  and  Esther  Ludington)  on  their 
farm  In  Idaho-  Her  married  Nancy  A-  Butterfield  in  August  1970.  After 
receiving  a B.S.  degree  in  electrical  engineering  from  Oregon  State 
University  in  June  197S  he  moved  to  Oainesville,  Elonida,  to  study 
astronoB^  in  the  graduate  program  of  the  University  of  Florida.  In 
NovesAer  1972  his  son  Hiilip  Jay  Ludington  was  born.  In  June  1978  he 
married  Karen  Faige  Rockwell;  three  years  after  a divorce  from  his  first 
wife.  Shortly  after  their  marriage  he  and  Karen  moved  to  Austin,  Texas, 
where  he  had  a postdoctoral  research  associateship  in  the  Departnent 
of  Astronomy  at  the  University  of  Texas  at  Austin,  The  degree  of 
Doctor  of  Philosophy  from  the  University  of  Florida  is  expected  to  be 
conferred  in  OeceiAer  1978. 
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